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1.1 Piecewise Smooth Functions A function f(x) is said to be piecewise : 
continuous in the interval [a, b] (the notation [a,b] will be used to ἀδηοῖθ 
the closed interval a < x < b) if there exist a finite number of points — 
X4,Xz,--+,Xp, (@ = X, < Χ... < Χ, = 5) such that f(x) is continuous in 
each sub-interval x, < x < x,,, amd has finite one-sided limits at the ἐν 
end points of each interval; these left- and right-hand limits will be denotes δ 
by f(x,+0) and f(x,,,—0) respectively. Clearly f need not necessarily — 
even be defined at the end points of the sub-intervals. 

A function f is said to be piecewise smooth in the interval [a,b] if fi 
and its first derivative f’ are both piecewise continuous in this interval. _ 


Problem 1.1 Show that, in the following, f(x) is piecewise smooth in its” ι 
interval of definition. Determine the points at which f or 7΄ are "εν 
continuous. 

(a) f(x) =x+e, («πες χε, f(x)=c-x, O< x <o). 

(Ὁ) fx)=0, (-1<x<0, 70) Ξ χ, O<x<t 


Solution. (a) f is obviously continuous in each of the intervals ve 
—c<x<0,0<x<c and f(0+0) = f(0—0) = c: thus f is con- 
tinuous for x in [—c,c]. For x <0, f’=1 and for x >0, f'= -1; 
hence f’ is continuous in each of the sub-intervals —c < x <0,0<x< aM 
but is discontinuous at x = 0. f’ is not defined at the origin but both — 
f'(0+0) and f’(0—0) exist and are finite (being + 1 respectively). Thus 
f' is piecewise continuous and f is piecewise smooth in [—c,c]. The Ὁ 
graph in Fig 1 represents ἢ v4 


a 


6) f is clearly continuous in each of the regions —1 < x < 0, 
0 <x <1and f(0+0), f(0—0) are both equal to zero. kis f’ = Ofor 


—1<x<Oand f’ = 1 for 0 < x < 1, so that {5 is discontinuous at 


x = 0. f’ is not defined at x = 0 but both f’(0—0) and f’(0+0) exist and 
are finite, being equal to 0 and 1 respectively. Thus / is a piecewise smooth 
function. oO 


1.2 Fourier Series If a function f(x) is integrable over the interval 


[a,a+2c] then the Fourier series of f on the interval is the trigonometric 
series 


1 = nx Μπχ 
= ττος In -— 1 
549+ δ, (4,00 ᾿ + b, sin ; ) (1.1) 
ἧς a+2ec 
where a, = =| f(x) cos dx, (1.2) 
a+ 2c 
b= : f(x) sin—— dx. (1.3) 


Fourier series are often defined over one of the intervals [0,27], [ --π, x], 
the results for these are obtained from (1.2), (1.3) with appropriate choices 


_ of aand c. The a, and b_ are known as the Fourier coefficients of ἢ, 


[An alternative form for the Fourier series can be obtained by expressing 
cos ntx/c, sin ππχίς, in terms of exp(+inzx/c) and this gives 


% + > a, exp(innx/c)+ Σ %, εχρί -- ἱππχ οὶ. 
πΞῚ a= 


_ where the bar denotes the complex conjugate and 


1 'a+2c 


a, = >. f (x) exp (—inzx/c) dx. 


This form can be further simplified by replacing n by —n in the second 
series, giving δ᾽ «a, exp(innx/c)]. 


8 838 «“)ῳ 


Theorem 1 If f(x) is piecewise smooth in the interval [a,a+2c] then 


_ the Fourier series of f on [a,a+2c] converges for all x in the interval. 
The sum, S(x), is equal to f(x) at all points where [7 is continuous; if / 


has a discontinuity at x = x, then 2S(x,) = f(x,+0)+ f(x,—0). At the 


end points 2S takes the value f(a+0)+ f(a+2c—0). 


Theorem 2 (Parseval’s theorem). If f is piecewise continuous in 
[a,a+2c] then 


’ ; - = "A ᾿ 7 a P| — a . = —~ . Ὗ <> , 
a .* ald ve i“ . ᾿ ᾿" 1 Ὁ νὰ a 
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1 fat 2c oo 
f? dx =4a+ Σ (a? - 83). 
Problem 1.2 (ἢ) Find the Fourier series of f(x) =x on [0,2π] αὖ 
determine the range of values of x for which the Fourier series converges 
to f(x). 

(ii) Sketch the function represented by the Fourier series in [ — 2π, 47]. 


Solution. (i) The coefficients a, and b, are found by setting a = 0 and 
¢ = 7 in (1.2) and (1.3). 


1 2x 1 2n 
a -=| x cos nx dx, b -τ| x sin nx dx, 
π π | 


1 2n 
= =| Χ ΟΧ = 22%. 


sae ἢ 2" sin nx dx, n # 0, 


ce 7 2 
x=2n 2a 
5, = ~| cos ns +— cos nx dx, 
x=0 0 
2 ὃ. 2 2 
=-——-+——sin2ine = —-. 
n nT n 


Hence setting c = 7 in (1.1) shows that the Fourier series of x in [0,27] is 
| = sin nx 
π-- 2 Σ ae 


n=1 
(ii) x is continuous and differentiable throughout the interval [0,27] 
and hence, by Theorem 1, the Fourier series converges to x at all interior 
points. The value of the sum at the end points of the interval is the mean of 
the limiting values of x, and so the sum at each end point is z. The Fourier — 
series has period 2z [i.e. S(x) = S(x+2z)] because each term has this 
period. Thus the sum S(x) represents a function of period 27, equal to 
xin 0 < x < 2xand to7z at x = 0, 2z. S is shown in Fig. 2, where dots 

denote values at x = —2z, 0, 27, 4π. 


Generally, if S(x) = f(x) for 0 < x < 2c, then, outside the interval, 
S(x) will represent that function of period 2c which is equal to /(x) on 
0 < x < 2c and the Fourier series can be regarded as representing, 


0 
Fig. 2 


outside the interval Ὁ < x < 2c, the periodic extension of f. The Fourier 
series cannot represent a function f(x) for all x unless f has period 2c. 


Problem 1.3 Find the Fourier series of cos 42x/c in [0, 2c]. 


Solution. In this case a in (1.2) and (1.3) is zero and hence 
1 | zc τῆν = 4x 
a=- cos —— cos —— dx 
c ς ς 


0 
2¢ = 
= =| cos = os FO | 
ΧΟ, ς 


1 [2 ππχ 42x 
b,=- sin —— cos —— dx 
a ci. ς ς 


edi | a 
= 


Forn # 4 the integrals are all zero and hence a, = b, = 0,n # 4. Setting 
_ n= 4givesa, = 1 and b, = 0, i, the Fourier series of f(x) = cos 4nx/c 
is cos 4nx/c. Oo 


> 


alt + on] ἊΣ 
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ἶ This result illustrates an important theorem: If two trigonometric series 


. -«- x 
Σ (a, cos +6, sin), Νὰ (« cos — +b, ,sin™=*), 
Γ n=0 . . n=0 : 
converge to the same sum at all but at most a finite number of points in 
[a,a+2c] then a, = a’, b, = δ΄. An equivalent statement is that if a 
Fourier series is identically zero then so are all the coefficients. 

Another corollary is that the Fourier series of any finite series of the form 


bh κ 
ΠΟ Σ [a, cos (nxx/c)+b, sin (nnx/c)] is the series itself. 
n=0 


Oo Be πε μόνον ts i I ΨΥ πόλου ee .: 


This latter result can be proved by direct calculation of the Fourier 
coefficients, the coefficient of cos mzx/c in the Fourier series is, by defini- 
tion. 


i ad _ mmx mmx 
- Σ a, cos +b , Sin —— cos —— dx. 


a=? 
The following results can be verified to be valid: 


Ζς. ἢπχ ΜΙΝΧ Ζ.ς. πχ =o MTX 
cos —— cos —— dx = sin —— cos —— dx 
ἃ ς ς P ς 


sin —— sin —— dx = 0, 
ὃ c c 


me sin? ex gx = 1 0, 

Ἀ ς 0,n = 0. 
Since the summation is over a finite number of terms the order of sum- 
mation and integration can be interchanged and the above expressions 
show that the only non-zero contribution is that for m = n giving that 


the Fourier coefficient of cos mzx/c is a,,(m < N) and zero otherwise. 
A direct consequence of the above analysis is that if f(x) is assumed to 


2¢ ππχ . mmx 
n#~ m, 


have an expansion of the form δ᾽ (a,cos(nmx/c)+b, sin(nzx/c)) then 
n=0 


carrying out the above analysis (formally replacing N by 00) shows that — ᾿ | 


a, and b, are the Fourier coefficients. This analysis is only formal as it his 
assumed a particular form for f(x) and also assumed that integration 
and summation can be interchanged in an infinite series. 


Problem 1.4 Find the Fourier series on [—1,1] of the function /(x): 
f(x) = 0, -1 < x < 0; f(x) = x,0 < x < 1. Hence show that 
Σ (»--1} 2 = 


r=1 


Solution. In this example a = —1 andc = 1. The function /(x) does not 
have the same form over the whole range in the integrals defining a, and 
δ... Thus the interval of integration has to be split into the two intervals 
[ —1,0] and [0, 1]: 


mf, 


— 
b= ‘a f(x)sin ππχ dx = ᾿ x sin ππχ dx. 


2 
I 


1 
f(x)cosnnmx dx = ᾿ x cos nmx dx, 


[f=0,-1<x <0] 


For n # Ὁ the integrals can be evaluated by integration by parts: the 


ἭΝ 


ΤΣ 


ἡ 


. i 


{ 


} 
y 


a 


case ἢ = 0 involves an elementary integral giving a, = 4. Thus 
᾿ : x=1 1 Ὡς 
ΠΣ Ed AST gin ges dy ΘΗΝ ΣΙ ete ἘΞ {|| Ὸ 


ΣῈ 
oe ΠΣ nm n>? 
που’ 1 [* 

| ————— +— | cosnnx dx 

| nn an So 

coom 1 . (—1)"*? 
 ——_ +; sin 2 = ---- «Ἐς 
nt 


Thus a, = 0 (n even), and a, , = —2/(2r — 1)?n?, r= 1, 2, .... The 


_ Fourier series is 


1 2  cos(2r—1)ax  (—1)'*'sinnnx 
4 π ᾿ς 51 ae > i 
Since f is a piecewise smooth function [Problem 1.1(b)] it follows from 
Theorem 1 that the sum of the Fourier series at x = 1 is 
al. f(I)+ f(—1)] = 5. Substituting x = 1 in the Fourier series gives 
Σ (2r—1)~? = 4n?. o 


r=1 


Problem 1.5 Find a Fourier series converging to |(sin x)| for every x. 


Solution. A Fourier series cannot represent a function for all x unless 


_ it has the same period as the function (in this case 2; see Problem 1.2). 
__ The precise interval is immaterial as both |(sin x)| and the relevant sines 
_and cosines have the same period, and the integral of a periodic function 


over a cycle is independent of the end points. This can be proved as 


follows: Given F(x) = F(x +2c) then (1-2), eee Hence 


[°° Fox) dx = |" F(x) dx+ |" *° F(x) dx = 


1] (* F(x) dx+ Γ΄ *° Ε(Χ 20) dx 
ὰ i F(x) dx+ | “ F(x) dx = i F(x) dx. 


Thus the integral is independent of a. In the present problem we choose 
the interval [0, 2] for simplicity, ie. a = 0, c = ἐπ in (1.2), (1.3). Hence 


a = = ὙΠΕΡ ΌΚΤΟΝ : 
ΝΡ τς χίδη -- 1) 


δ = = | sin x sin 2nx dx = 0. 


By Theorem 1 the Fourier series of |sinx| converges to |sin x| at all 
6 


wins % “34 ΡΨ a J ; δὰ ἰδ. Ὁ 
“a. 7 @ 7 y-e-7 ἀξ ν᾿ τος 


interior points of [0,2]. At the end points the series converges to the 
mean of the end values of |sin x|, i.e. to zero. It then follows by the periodi- 
city of |sin x| that, for all x, 


᾿ 2 4 5 — iH 
jsin x] = ---- = Σ τὰ ᾿ Oo 


Problem1.6 If f(x) = sinx,0 < x < x, f(x) = 0,2 < x < 22,findthe 
Fourier series of f(x) on [0, 27]. 


Solution. This problem can be solved most easily by noticing that 
f(x) = {sinx+|sinx|). The Fourier series now follows directly from 
Problem 1.5 and the fact that sin x is its own Fourier series. The required 
Fourier series is 


j ὦ ~ cos 2nx 
Ἅ οἱ = 
aaa om π 2, (4n?—1) 0 


Problem 1.7 Find a Fourier series converging to sin? 4x for all x. 


Solution. sin? 4x is a function of period 2x and, as in Problem 1.5, it 
suffices to determine its Fourier representation in [0,27]. This representa- 
tion can be obtained most simply by writing sin?4x as 4(1—cos x) and 
noticing that, by the theorem quoted in Problem 1.3, the Fourier series 
of this function is itself. oO 


1.3 Fourier Cosine and Sine Series In certain circumstances the Fourier 
series of a function f(x) in [—c,c] consists only of sine terms or only of 
cosine terms. This arises when 


(i) f(x) is an even function of x. [i.e. f(x) = f(—x)]—cosine series 
(ii) f(x) is an odd function of x. [i.e. f(x) = —f(—x)]}—sine series. 
This is illustrated in the following problem. 


Problem 1.8 (i) Show that the Fourier series on [—c,c] of an even 
function f(x) contains no sine terms. (ii) Determine the general form of 
the Fourier series when f(x) is an odd function. 


Solution. (i) Setting a = —c in (1.3) gives 


cb, = | f(x) sin— dx. 


The range of integration is now split into two ranges [ —c,0] and [0,c], 
7 


" ᾿ > αι Δ Μ 
Ὁ» fa } sy 
my al «ὦ 
a "5" ᾿ 
ima a? © 


᾿ and x is replaced by —x in the former. Thus 
fo c 
ἥν ον, =| f(x)sin —= dx+ | f(x) sin dx 
dy | 


+e δ δ 
¥ = -| f(—x)sin χα] f(x)sin~ dx = 0, 
7 ἃ ς = c 
since f(x) = f(—x). Therefore the Fourier series will not contain any 
sine terms. From (1.1) we note also that 
0 ς 
ca, = | f(x) cos = 4.4} f(x) cos == dx 
πο" 0 


-| f(—x) cos — dx+ | f(x) cos — dx = 2 f(x) cos —— dx. 
, 0 c 0 c 0 c 


we (ii) For the case when f is odd, it follows from the above expressions 
a for a, and b,, by replacing f(—x) by — f(x) that a, = 0, 
- cb, = 2 f§, f(x)sin(nnx/c) dx, whence the result follows. " 


Ρνοδίοηι 1.9 Show that, in [ -- π, x], 


— (-1 
x? = tn? +4 » am 
n=1 
oe οο [4] 
and hence evaluate (i) δ᾽ π΄, (ii) Σ᾽ (—1)"*!n~?, (iii) ¥ n+. 
n=1 n=1 a=1 

Solution. The first step is to determine the Fourier series of x? on 
[--π.π]. Since x? is even the series will consist only of cosine terms, and 
setting c = x in the appropriate equation of Problem 1.8 shows that the 
Fourier series is 


ἵν 9] 2 ὄπ 
βαρ Σ᾽ a.cosmx where a -- | x? cos nx dx. 
π 
0 


n=1 


The integral can be evaluated by integration by parts for n # 0, giving 

n’a, = 4(—1)"; the integration for n = 0 is trivial and gives 3a, = 2n?. 

Since x? satisfies all the conditions of Theorem 1 and takes the same value 

at the end points of the interval [ —7, 2], the sum of the Fourier series is 
equal to x? at all points in [—7, π], ice. 


x? = 4n7+4 Σ᾽ [(-1)'cosnx/n’], [x] « π. 
n=1 


(i) Setting x = 7 gives, since cosnz = (—1)", δ᾽ πη = tn’. 


n=t1 


= 


(ii) Similarly setting x = 0 gives Py (—1f*4 n-? = ta’. 


οῦ 


(iii) The determination of Σ᾽ π΄ is less straightforward, but canbe 


n=1 


made by noticing that the terms are squares of the coefficients of the 


Fourier series for x”. This suggests use of Theorem 2 (Parseval’s theorem) 


which gives 
1 [5 3 πὸ ες 
-- x* dx = —+16 ee 
a ὁ 18 neh 


2 


Hence 900 Ὁ nt = x. Gta 


n=1 

For any function f(x) defined in [0,c] one can construct an even 
function h(x) equal to f(x) for 0< x <c by setting h(x) = {{-- χ) for 
—c < x <0 and thus, by Problem 1.8, there will exist a Fourier series 
representation of f(x) in 0 < x <c consisting only of cosine terms. 
Similarly by defining h(x) by — f(—x) for —c < x < 0 it will be possible 
to construct a Fourier series representation of f(x)inO0 < x < cconsisting 
only of sines. These representations are formally defined as follows: 


Fourier cosine series If f(x) is integrable on [0,c] then the Fourier cosine — 
series of f on [0,c] is defined to be 


Ὃ 
=a,+ ) ie 
2° ᾿ Cc 


n=1 


where a, = : f(x) cos dx. (1.4) 
0 


Fourier sine series If f(x) is integrable on [0,c] then the Fourier sine 
series of f on [0,c] is defined to be 


x b sin aot, 
ς 
n=1 
Σ fy 
where b == | f(x) sin — dx. (1.5) 
C's ς 
Theorem 3 The Fourier cosine and sine series of a function f(x) which 
is piecewise smooth in [0,c] converge to f(x) at any pointxin0< x< ¢ 
where f is continuous. At any interior point x = x,, where / is discontin- 
uous, both series converge to 4[ f(x, +0)+ f(x,—0)]. The cosine series 
converges to f(0+0) at x = 0, to f(c—0) at x = c and is always con- 
tinuous. The sine series converges to zero for x = 0, c. Theorem 3 can be 
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proved i 

cs equal to f(x) in 0 < x < c. (Since these series only represent 
on an interval whose length is half the period of the functions involved, 
the series are often referred to as half range series.) 
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by applying the results of Problem 1.8 to the even and odd wie 
x). 
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τ (b) The coefficient b, in the sine series of cos }x is 


i? es Se | 
b= : sin nx cos$x dx = ἊΝ [sin(n -Ὁ Σ)χ + sin(n—4)x] dx 


7 " Ν 
pi Tt should be noted that a corollary of Theorem 3 is that for a function 5 εὖ Se | ‘ 
~ f(x) continuous in [0,c], the Fourier cosine series converges ἰ ΓΝ aes. me 
this interval but the convergence of the sine series in the closed garb The series converges to cos 4x at all interior points of [0,7] and to Ζεῖο 
‘requires the ad ditional conditions [(0) = f(e) - 0. Another ἔρνη ag! at the end points and hence represents cos$x in0 < x < x. ‘ 
“ is that a function continuous In [0, ¢] and vanishing at x = 0 and x = ς The weriendelines an‘ οὐδ fanction of x and hence in sao) eae ‘4 
can be represented throughout the interval by its Fourier sine series. oan μι costx. From this and the periodicity property it follows that 
“hs : ᾿: the graph is as shown in Fig. 4 (the form of the series in [ — being 
Problem 1.10 Find, on the interval [0, x], (a) the Fourier cosine seri π,π] being | 
of x and (b) the Fourier sine series of cos 1x, Sketch the functions defined denoted by heavy lines). a sa 
ἷς ΒΥ the Fourier series for |x| < 3m. ἤ 
. . Solution. (a) The Fourier cosine series of x is 5a)+ py a, cos nx where “ἂ 
Εν π -- i ion by parts - 
᾿ς αᾳ = (2/n) [5 xcosnx dx. We find a, = 7, for n # 0 integration “ 
oa τᾷ Jo 2[(—1)"—1]/x. Thus a, will vanish for even values of n. For 
sre) (= 1,2)...) we get x@r—1)'a,, = 4. The Bours | 
cosine series is thus “a 
ΧΆ 4 — cos(2r—1)x ὯΝ 
— ἘΠΕ (2Ζγ- 1} ie 
and by Theorem 3 converges to x in [0, 7]. . é Ὁ 
ο΄ ΒΕ cosine series is an even function of x and Is equal to x in [0, x] and to a 
Ra cin [ — n, 0]. This even continuation (denoted by heavy lines) of x is a 
a shownin Fig. 3 together with its periodicextension (period 2n)to [ — 3n,3n]. re 
ἢ f 4 | ἥ ᾿ 
Α ὯΝ Fig. 4 ad 
Ta A! 
Via Problem 1.11 (i) Determine the general form of the Fourier cosine and 
ee π τῶ series in [0,c] of a function f(x) satisfying the relationship f(c—x) = 
. αὶ, 
| | (ii) Obtain the corresponding series when f(c—x) = —f(x). ‘a 
Solution. (i) We consider first the integral from x = 0 to x = c of an | p 
; arbitrary function F(x) such that F(c—x) = F(x). Such a function is 
. : =< ama κ᾿ ae τ το said to be even about 4c as setting x = 4c—« gives Féc+e) = Fbc—e), 


ce i.e. it takes the same value at a given distance on either side of x = 4c. ia 


A Fig. 3 [If F(e—x) = —F(x) then the function is said to be odd about x = 4c.] Ἷ d 
90 “͵ 
ἜΝ Β ν᾽ 
εν... ag 
ἯΙ “4 
Pr e τ 
Li Ἂν ‘ ὶ 


᾿ 


Ν᾿ 


e ͵ 
"ἃ 
1 
τὰ, 5 
7 ε"ς _ ‘ 
' 
“Ἂ r ᾿ ᾿ ι 
ἐ Ν i ye a 
Z 


_ The range of integration is split into the two intervals [0, 4c], [4c, c], thus 


4 ἃ a -* 
2 ᾿ 
¢ s 
5 Ἵ 


[ F(x) dx = ᾿Ξ F(x) dx+ [: “ F()dx 


= |" F(x) dx+ ey F(c—x) dx = 2 [ἢ FG) dx. 


For a function G(x) odd about x = 4c it follows that [6 G(x) dx = 0. 
_ Sinnmx/c is odd or even about 4c according as to whether n is even or 


odd and cos nmx/c is odd or even about 4c for according as to whether n 
is odd or even. Thus, for the given fonction f(x), f(x)cos nrzx/c is even or 


" 


a 
{ 


odd about +c according as to whether n is even or odd and from the above 
general result it follows that a, = 0 in the cosine series for n odd and 


4 (te 
a, Ἐπ᾿ f(x) cos ax, y= O:1,.... 


be ἊΝ (χ) 5ἰπ ππχίς is odd or even about ἔς according as to whether n is even or 


_ odd and thus b, = 0 for n even and 


b,,-; -: Πρ) sin (2r—1)— ἀν, Pm LBs 


(i) Im this case f(x)cosnnx/c is even or odd according as to whether 


_ nis odd or even and thus a, = 0 for n even and 


Gs - f(x) cos (2r— n= dx, fm 1.22... 
"Similarly b, = 0 for n odd and 


a τ" f(x )sin™ dr, Yee oO 


᾿ Problem 1.12 Find the sine series for cos x in [0,π]. 


a” 


Solution. Cos x is odd about ἐπ and hence, from Problem 1.11, the sine 


series will be of the form >) b,, sin 2rx where 
r=1 


8r 
m(4r2 —1) 


The sine series converges to zero at x = 0 and x = mand hence 


4 (i 
b,, = = | cos x sin 2rx dx = 
eT) 0 


coxa ty ἢ τ: ΑΙ O<x « π. oO 
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Replacing x by 1x gives A ‘i 
8 τ rsinrx ᾿ς ΡΝ 
iL ee 7 ‘< 4 y 
C08 3x = — a Τ᾽ O0<x ΖΗ, a. 


This latter result is in Jat with one previously obtained in Problem — 
1.10. 


Problem 1.13 Find the sine and cosine series on [0,c] of the function — 


f(x): f(x) = x, 0 < x < 4c; f(x) = c—x, ἔς < x < cand hence evaluate 


a ᾿ 


Σ᾽ (2r—1)~?. Sketch the functions defined by the Fourier series in 


ix! < 2. Ἂν 


ν᾿ Ὃ 
Solution. The cosine series ἰ8 Σὰ. ἘΣ a, cos(nnx/c) where ae 


n=1 


fa nmx 
= =| J (x) cos —— dx. 


" 
᾿ 


a 
- 


ἱ 


f(x) does not have the same form over the whole interval of integration _ 


and this interval has to be split into the two intervals [0,4c], [4e,c]. fis — 
even about x = $c and from Problem 1.11 it follows that a, = 0 for odd 
values of n. However in order to illustrate the point that the correct 
Fourier series can be calculated without direct use of the evenness property — 
the integral for a, will be —er directly. 


ca te 
= | x cos —— dx + " (e=x)e0s = dr, n#0 


2 
- 3 2 cos — | —(- ir | n # Ὁ (by integration by parts), 


ay = te. 
Clearly a, vanishes for n odd and hence 


c 
a,,= πΣ,2 (cos rm— 1). 


an 


aa 


ΔῊ 


0 ἐς ΠΣ 


This vanishes for r even and the only non-zero values are forr = 2s—-1 ὁ ὁ ὁὃϑΘὃὅ Γ 


[s = 1,2,...] when a,._, = —2c/n?(2s—1)?, the Fourier series is thus 
1. 2ς “ cos(4s—2)nx/c 

4΄ πἪ oS, ὦ- 1) 

f(x) is piecewise smooth throughout [0,c]. Hence, by Theorem 3, the 

Fourier series converges to f on [0,c] and setting x = 0 gives 


y (2s—1)~? = 1π2. 
13 


Ji ἡ ie a Ὧ 


οἱ ΕῚ 


—2c =e 0 Cc 2c 


Fig. 5 


The cosine series represents an even function of x equal to f(x) in [0, c] 
and its form is shown in Fig. 5. For c = 2z this has the same form as the 
Fourier series of Problem 1.10 (a). 

The coefficients b, can also be calculated directly as above or by using 
the evenness property of f(x). Problem 1.11 shows that b,, = 0 (n even) and 


$c 
cb, _, = 4 x sin(2r—1) ax. r= 1,2,.... 
r ὃ Cc 


— 4-19" ᾽ς 
ΝΠ ὯΝ 
f(x) vanishes at x = 0 and x = c and thus, in [0,c], 


[by integration by parts]. 


ἄς oo (-1)"' . 1X 
f(x) = a) ‘a ἔπεσ: ὩΣ De 


The sine series represents an odd function of x equal to f(x) on [0,c] and 
of period 2c: its form is shown in Fig. 6. 


r=1 


Fig. 6 


ΒΥ ΡΥ. 7 . δ a 
’ i‘ ᾿ : * .. 

- 4 : 5 -" ν ) ᾿ ᾿ ᾿ἷ > ion J " : Ι = i “ 5 ᾿ 4 
᾿ τὰ - * ' ie ᾿ * Me, a ἢ 


1.4 Cosine and Sine Series of Period 4c In applications it is often 
necessary to represent a function f(x) on [0,c] as a cosine or sine series ὖ 
whose general terms are a,,_ , cos(2n—1)mx/2c and b,,_ , sin(2n—1)nx/2c 


respectively. From Problem 1.11 it is seen that such representations could 


be obtained on [0,2c] for functions odd or even about x = c, and thus all © 
that it is necessary to do is to extend the given function to [0, 2c] and make 
it even or odd about x = c as necessary. The coefficients a,,_, ,b,,,_, can 
be obtained by replacing c by 2c in Problem 1.11, giving 


7a kg 2n— 1)πχ 
a2,-1 = al ΓΟ ρου “5: ΤΣ dx, 


2n—1)nx 


2 ie 
Be = τ | J (x)sin Φ-ΞΞ ax 


The even extension of f is continuous and thus the sine series converges 
to finO < x < candtozeroat x = 0; the odd extension is discontinuous 
at x = c and the cosine series converges to f in 0 < x < c and to zero 
atx ΞΞ 6. 


Problem 1.14 Determine the Fourier cosine series and sine series of 
period 4z which represent x in [0,2] and sketch the functions represented 
by the series in |x| < 4z. 


Solution. Setting c = x in the preceding analysis gives 


@ 
x= >} a,,_, coss2n—l)x, O< x «κα; 
n=1 


x= y b,,,-, sin3(2n—1)x, 0<x<Q, 


where 
ΣΡ 2 a 
a.,-1= =| x cos>(2n—1)x dx, b,,_, = al x sin4(2n—1)x dx. 


The integrals can be very easily evaluated by integration by parts giving 
Mayet 8 8(-1}7 
τι 0η-ἢ πρπ-Ὄ . πὶ ΝΡ. 

x vanishes at x = 0 and hence the sine series converges to x in [0,7]. 
Cos 4(2n—1)(n—x) = —cos4(2n—1)(n+ x) and thus the cosine series 
represents a function odd about x = 2 (this of course also follows from the 
discussion preceding this example). Thus the form of the cosine series is 
known for 0 < x < 2x and, as the cosine series is an even function of 
x, it can be now determined for --- 2ῶπ < x < 0. The series is thus known 
for a complete cycle and the form for all x can now be obtained and is 
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1.5 Differentiation and Integration of Fourier Series Differentiation of 
a Fourier series will produce another Fourier series whose coefficients are | 
proportional to na,, nb, where a, and b, are the coefficients of the original 
series. Thus, even if the original series converges, the differentiated series 
may not even converge if a, and b, are not sufficiently small for large n 
and it is therefore useful to be able to derive, from the form of the original 
function, the behaviour of the coefficients for large n and we have the 
following theorems: 


| Theorem 4 If f is piecewise smooth in [a,a+2c] then a, and ῥ are at 
Fig. 7 : most O(n~*) (ie. lim na,, lim nb, are both finite). 
ey ; ΔῈ , Theorem 5 If f and its first (p—1) derivatives are continuous in 
given in Fig. 7. The form for the sine series is obtained similarly and is [a,at+2c]; fa) = f(a+20), r = 0,1,...p—1, (where f is the rth 
given in Fig. 8. " derivative οἵ f) and f is piecewise smooth in [a,a+2c] then a, and 
b are at most O(n~?~'). 


Problem 1.16 Show that, for large n, the coefficients of the cosine series 
on [0,2] of a function f(x), continuous with a piecewise smooth first 
derivative, are O(n™ ”). 


Solution. The behaviour of the coefficients cannot be deduced directly 
from the above theorems which refer to the complete Fourier expansion 
of a function. The cosine series expansion of f is the complete Fourier 
expansion on [—7,7z] of the extended function even in x and equal to 
f in [0,2]. This function is continuous in  -- π, 2], taking the same value 
at both end points, and hence Theorem 5 may be applied with p = 1, 
Fig. 8 : giving that the Fourier coefficients are O(n™ ?). oO 


The derivative of the extended function is not continuous at x = 0 as 
These results provide two more examples of Fourier series expansions | its left and right hand limits are + (0) respectively, at the end points this 


which represent the same function over a given interval but which repre- derivative is +f'(m). Thus, if f(0) = f'(x) = 0, Theorem 5 show? that 
sent very different functions when considered for all x (cf. Problem 1.13). the Fourier coefficients of the cosine series for f would be O(n~°) if f 
were continuous and f{” piecewise smooth. 


Problem 1.15 Find the representation in [0,c] of sin zx/c as a Fourier Problems 1.10(a), 1.13, 1.15 are all examples on deriving cosine series 
cosine series of period 4c. of functions satisfying the conditions of the present example and all have 

Solution. The coefficient a,,_, of cos (2n—1)nx/c in the expansion is Fourier coefficients of O(n-*). It follows that the coefficients of the sine 
given by series obtained by differentiating cosine series of the above function are 


O(n~*) and the series is at least convergent. This however does not prove 
᾿ = gh τ τὸς πχ cost Qn— Ha ae that the differentiated series converges to f(x). Similarly the coefficients 
2n-1 : 2 ers of the Fourier sine series of f are O(n~*) but will be O(n~ *) provided that 


The integral can be evaluated by expressing the integrand as a sum (0) = f(x) = 0. Thus in many cases the series obtained on differentiating 
of sines and this gives a. = —8/nx(2n+1)(2n—3). 0 | the sine series of a continuous function will not even converge. 
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The general conditions concerning differentiation of a Fourier series 
are contained in the following theorem: 


Theorem6 If f(x) is continuous and f(x) piecewise smooth in [0,c] then 
(i) The Fourier cosine series for f may be differentiated term by term 
and the resulting series converges to f’in0 < x <<. 
(ii) The Fourier sine series may be differentiated term by term if, and 
only if, (0) = f(c) = 0, and the differentiated series then converges 
to f’ in [0,c]. 
For the case when /{(0) or f(c) are not zero the coefficients of the cosine 
series for ζ΄} are given by 


a= =[f0-fO], 4, = 


In applying this theorem to cases be ΤΣ derivative of f is not defined 
at a point x = x,, f’ at that point is to be interpreted as 
SL f'(x,+0)+f'(x,—0)]. The theorem is proved by using integration by 
parts to express the Fourier coefficients of f’ in terms of those for ἢ 


Problem 1.17 Use the results of Problem 1.10(b) to obtain the Fourier 
cosine series of βίη ἔχ on [0, z]. 
Solution. Problem 1.10(b) shows that 


nnb, 


mal 8 = nsinnx ὁ «νὸς 

--ὀπὁἰς == pale x <7. 

2 π pat mde, ἣν 
n= 


The Fourier coefficients a, of —}sin 4x are, from Theorem 6(ii), 


᾽ “»᾿ , 8η2 eet 2 
a og n 7(4n?—1) π 2(4n?—1)’ 
2 4 Ξὶ cosnx 
Th Ae pc eng, ne ὦ, ὡμ , 
us sin3x = — ἘΣ ato) O<x<n O 


Problem 1.18 Use the results of Problem 1.13 to determine the cosine and 
sine series in [0,c] of the function g(x): g(x) = 1,0 < x < 4c: gc) = 0; 
Ax) = -lic<x<€e. 


Solution. Except at x = ἔς the function g is clearly the derivative of 
f defined in Problem 1.13. At x = 4c, however, f’ is not defined and thus 
7΄ at x = 4c, for application of Theorem 6, has to be interpreted as the 
mean of the limiting values of the derivative on either side of x = 4c, 
ie. zero. Theorem 6(i) then shows that the sine series for g is 


ay [sin(4s — 2)nx/c]/(2s — 1). Theorem 6(ii) shows that the coefficients 
s=1 
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in the cosine series of f’ are given bya, = O[ f(0) = f(c) = 0], a, = nnb,/c, 
where 4(— γ 1 


b,=0 (n even), b5,-1= Qr—1) m2 r= τον 


This gives a, = 0(n even) anda,,_, = 4{-- 17 */2r—1)a,r = 1,2.... Π] 


Results analogous to those of Theorem 6 can be obtained for the cosine 
and sine series of period 4c defined in section 1.3. 
Theorem 7 If f is continuous in [0,c] with a continuous first derivative 
then 
(i) when f is represented in [0,c] by the series 
= πχ 
& a, , cos(2r— ) τς 


the derivative /' can be represented by the series 


> πχ 
rs ba, 1 Sin(2r—1) το 


+1 
Birks ee : Han στ 


(ii) when f is represented i in [0,c] by the series 


πᾶ;͵ .- 1° 


», b,,-, sin (2r— > 
the derivative f’ can be represented by the series 
Σ ἂλ,: αὶ cos(2r— ἢ == 
(2r—1)nb,,_, 
2c 
This theorem is also proved by means of integration by parts. 


2 
where a,,-4 ΤΊ TIS 


Integration of Fourier series The situation concerning the integration 
of Fourier series is much simpler than that for differentiation and it can 
be proved that a Fourier series of a function f may be integrated term by 
term and that the integrated series converges to the integral of f. 


Problem 1.19 Show that the Fourier cosine series of sin }x on [0,2] may 
be obtained by integration from Problem 1.10. 


Solution. In Problem 1.10 it is shown that 


costt By = sin nt O<t<n 
: Pa or a 


“ΘΝ win iyo = 0 tot = xgives_ ἈΝ ΣΡ | pie Sin ἐπέ ἮΝ Fewer cosine 
| - nx -- ἸῈ a ἔν of (cr function g(x) defined in Problem 1. 18. 
na = --ς᾽.. ——— ἀπε; ΓΞ: ἀξ: aah” 
ΣΝ Dee -1 ’ sa γ ἧς: 4. Find the cosine and sine series of period 4c representing exp a Ἂν 
0<x<c. 


Si _ 5, Determine the cosine series in [0,2] for sin x by applying Theorem 6 δ᾽ 
cos nx | to the sine series for cos x derived in Problem 1.12. 


An’ 1 


᾿ I It is known that Σ (4n?—1)~' = 1 and thus 


| The sum Σ (4n?—1)~? saith lig ee = 0 in the co- 
- sine Mielec for sin xin [0, x]: 
i sinx = 7-= ) pa ical O<x<n O 


. an =I 
<< 
gil _# Ὁ the integrated series + wil not be a Fourier series. 


EXERCISES 
1. Show that 


ae 
. @ <x < 2. 


a=1 n=1 
2 sinh az ταν (—1)"acos nx —nsin nx 
π 2a a? +n? i 


n=1 


(ii) exp ax = 
4 


—h<X<7. 
By considering the value of the series at the end points show that 


1 «= @ 
jnacoth πα = ΞῈ Dany 


ws n—x)Qn—x)= δ᾽ =F O<x ς 2π. 

n=1 ὼ 
cos(2n — 1)x 

(2n—1)* ’ 


τ ra 
ere) 


n=1 


XS, 


and hence evaluate δ᾽ (2η-- 1)" 2, Σ᾽ (2n—1)~¢. 
n=1 n=1 
__ Determine on the interval [0,7] the sine series for x and the ‘cosine’ 
_ serie es for cos}x and obtain the range of values of x for which the series 
e “μά ἴο the relevant functions. 


. 


Chapter 2 


Solution of Boundary-Value Problems by Means 
of Fourier Series 
2.1 Introduction Most of the subsequent problems will be concerned 


with the solution of the simpler partial differential equations of mathemat- 
ical physics. Particular examples of these equations are: 


V7u = p (Poisson's equation), 
1 67u 
V7u = ot at (wave equation), 


V7u = ku (Helmholtz’ equation). 


V*u = 0 (Laplace's equation), 


V7u = Kk (diffusion equation), 


For two independent variables (which is the case with which we shall 
be concerned) all the above are particular cases of 


ας δα ον ᾿ς Ἐ ἀμί) + by Oy oa +b) tbs = ΤῸ») (2.1) 


known function of these variables. 

The methods to be described are only suitable for solving (2.1) in a 
rectangular domain D of the x, y plane and in the present chapter it is 
necessary that the range of one independent variable is finite and that the 
coefficient of the second derivative with respect to that variable does not 
vanish within D. Without loss of generality x is taken to be that variable 
(i.e. a, #0 in D) and D to be the domain 0 < x < Ἢ b < y < d; either 
or both of b, d can be unbounded. 

The boundary conditions under which (2.1) is soluble depend critically 
on the sign of b, /a, (it will be assumed this sign does not change within ἢ). 
[An equivalent method of relating the boundary conditions to the 
equation is by reference to the classification of partial differential equations 
in two independent variables (cf. P. M. Morse and H. Feshbach, Methods 
of Theoretical Physics, McGraw-Hill, 1953, Ch. 6). In this nomenclature, 
___ b,/a, being positive, negative or zero determines whether (2.1) is elliptic, 
᾿ς hyperbolic or parabolic. ] 

__ The boundary value problems to be solved in subsequent problems 
will be such that they possess unique solutions and the following table will, 
for most practical situations, indicate the appropriate conditions to be 
applied. 
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a Se τὼ Te? a. >) ie Te. oe Ἢ ἢ 


b,/a,,+ve πογ δι ον givenony=b,y=d 
b,/a,,—ve μ, δι ὃν given on ν = b. 
b,/a,,9 u given on y = ἢ. 


gk ao Cu 23 u Ou 


where a,, b, are known functions of x and y respectively and f is another - 


ee 


u or Gu/éx given 
for x= 0;c. 


2.2 Homogeneous Problems Homogeneous problems in the present 
context are defined to be those of solving 


a, a7 +4, +4 yu+b, 5 Ba ag, Oe (2.2) 
in D with u or du/dx vanishing on x = 0, c. The use of Fourier series (and 
their generalizations) to solve (2.2) relies on two separate concepts. The 
first is the principle of superposition which states that if any number of 
functions separately satisfy (2.2) then so will any linear combination of these 
functions (this is because the left hand side is linear in u). The second con- 
cept is that of separation of variables which consists of seeking product 
solutions of (2.2) of the form X(x)Y(y). The basic technique in subsequent 
problems is essentially the determination of combinations of product 
solutions which are such that they satisfy all the conditions imposed. 
This approach is generally known as Fourier’s method and is illustrated 
in the succeeding problems. 


Problem 2.1 Solve the equation é7u/éx?+07u/éy? = 0, 0 «χ « ς, 
0 < y < d,under the boundary conditions u(0. νὴ = πίε, vy) = 0,0< γ <d, 
μίχ. 0) = f(x), u(x,d) = 0, 0 <x <c, where (i) f(x) = βίη πχϑο, (1) 
f(x) = x(c— x). 

(One physical interpretation of this problem would be the determination 
of the steady state temperature in a rectangular region where three sides 
of the rectangle are kept at zero temperature whilst the fourth is main- 
tained at a known temperature f(x).) 


Solution. Assume u is of the form X(x)¥(y) where X and Y are functions 
of x and y respectively. Substitution of this in the equation for u gives 
a » 
Fi es ὅς 
the dashes denoting derivatives with respect to the relevant variable. A 
function of x can only be equal to a function of y if both are constant. 
Hence X” = —AX, Y” = iY. where A is a constant (the separation con- 
stant). The condition u(0, y) = u(c, y) = 0 gives X(0) = X(c) = 0 and the 
general solution for X is Ax +B, 4 = 0; Ccos/*x+Dsin A*x, 2 # 0. For 
λ = 0 the conditions X(0) = X(c) = 0 show that A = B = Ο and, for 
λ τὰ 0,X(0) = OgivesC = Oand X(c) = Othen requires that D sin λῆς = 0. 
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ie i od 


4 - ΗΑ , ν t 


ἂν», : 
4 2 “a 


ε΄ ee ee 
‘hus in order to have non-trivial solutions (1.6. D # Oca must be equal 
ἘΠῚ: nz, where n can be any integer. X will then be proportional to 
sin ππχίς and the equation for Y can be solved showing that Y will be an 
| - arbitrary linear combination of coshnzy/c, sinhnzy/c. Thus possible 
τ΄ solutions for u, satisfying the conditions at x = 0 and x = c¢, are 
Pia: cosh nzy/c+B, sinh nzy/c]sin nzx/c where n is any integer and the 
_ arbitrary constants A,, B, may depend on n. 
ο΄ In accordance with "the principle of superposition we thus consider the 
general form 
* n= > 4, cosh "+B, sinh— ml, sin 
ΧῈ n=1 
and attempt to find A, and B_ to satisfy the conditions at y = 0 and 
ΟΡ = d. The condition on y=d ‘will be satisfied if 


NIX 


A, cosh = = — B sinh = 


J and, for case (i), the condition on y = 0 gives 
a sin— = Σ᾽ 4, sin—, θ εχ ἃ 
C7 ot ς 

Ὁ Thus the A_ are the coefficients of the Fourier sine series in [0, c] of sin xx/c; 
__ this latter series is the function itself [cf. Problem 1.3] and hence A, = 1, 
A, ="0,n 4 land u = — cosech (zd/c) sinh πίν — d)/c sin πχίς. It ae be 
- verified directly that this satisfies all the conditions imposed. 

The only difference in case (ii) is that A, will now be the coefficients 
_ of the sine series representing x(c—.x) in [0, ch, i.e. 


sd 


q 2{° , 4c? 

= | A, -2 x(c—x) sin dx = 3,3 L1—cos nm], 

᾿ς Onintegration by parts. Thus A, = O(neven)and A,__, = [8c?/x°(2r—1)°], 
_ r= 1,2,.... The complete solution is thus 


. 8c? = 1 (2r—1)nd . , (2r—1)x(y—d) 
| 5 rs -- -τ pa pea ee ee eee 
Με (2r— Dex 


The fact that each term of an infinite series satisfies a differential 

' equation does not necessarily mean that the sum satisfies this equation, 
as this would imply that differentiation of the series produces another 
_ convergent series. This is not always the case [cf. remarks following 
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a ὧν ἐν κὸν 
Tee 
priiens " 16]. In the above series, however, the coefficients tend to Ἦν 


» «a 


zero sufficiently quickly for the series to be differentiated twice and the, “. 


sum does satisfy Laplace’ s equation. The boundary conditions are all + 


solution for case (ii). oO πω 


Problem 2.2 Solve the previous boundary value problem when — 
f(x) = x(c—x)and the condition on x = Ois replaced by u(0,y) = y(d—y), 


O<y<d. ἢ 


Solution. The boundary conditions on neither x = constant nor y= 
constant are homogeneous but, as Laplace’s equation is linear, u can 
be written as u,+u,, where u, and u, both satisfy Laplace's equation 
and satisfy homogeneous conditions i in one or other of the variables. We — 


have ᾿Ξ 


u,0,y) =u,(.y)=0, O<y<d, 


u,(x,0) = x(c—x), u,(x,d) = 0, O< x <4 
u,(x,0) = u,(x,d) = 0, - ae =e; . 
u,(0, y) = γί -- γ), u(y) =0, O<y<d. 


u, is the solution in Problem 2.1, case ἮΝ The problem for u, is very — 
similar to that for u, except that homogeneous boundary conditions are 
prescribed on y = 0 and y = d. A solution could be obtained by a very 
similar procedure to that of Problem 2.1 but in this particular case all _ 
additional work can be avoided by noticing that interchanging x and y 
and c and d reduces the problem for u, to that for u, . Thus the solution 


can be obtained from the previous result as a sum of two infinite series. 


This technique of splitting the problem into two separate ones with 
homogeneous conditions to be satisfied on one family of parallel lines 
in each case may be employed for any linear elliptic equation when 
non-homogencous conditions are prescribed on the boundary of ἃ ‘ 


rectangular region. ΠΚΠ on 


Problem 2.3 Solve the equation 


eu 1 Ou δι 

a pat ee Foe Se ed .-_ 

jee 0? OF" “2 
where v is a constant, under the conditions a 


u(0,t) = 0, du{x,t)/ox =0, x=1; t>O9, 
u(0,x) = x, éu(x,)/ét =0, t =0; O<x<l. 


(This problem can be interpreted physically as determining the longi- 
tudinal vibrations of a rod of length / clamped at x = 0 and free ai x = |. 
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Initially the longitudinal displacement at a distance x from the fixed 
end is equal to x and the rod is at rest. These initial conditions could be 
produced by applying a force at the free end parallel to the rod and remov- 
ing this force at time t = 0.) 


Solution. Assuming a product solution X(x)T(t) and proceeding as in 
Problem 2.1 gives Χ'' = —AX, T” = —Av?T. The conditions on X are 
X(0) = Χ ἢ = 0 and the general form of X will be that of Problem 2.1. 
For A= 0 the conditions X(0) = Χ΄ (ἢ) =0 again give A = B=0; 
for 2 τέ 0, X(0) = O gives C = 0 whilst Χ (ἢ = 0 requires that D cos J#/ = 0. 
Thus for a non-trivial solution 4*/ = (2n—1)n/2, n = 1,2,..., where n is 
any integer and the corresponding X is proportional to sin (2η -- [)πχ. 2]. 
Hence solving the equation for T and using the principle of superposition 
suggests the general form 


4, cos ata 1)xvt ᾿ ἷ Cnet με (2η -- [πα 


- Β 
u sin 5 


2] a, 2] 
n=1 
The condition δι δὲ = 0, t = 0, is satisfied by taking B, = 0 and 
u(x,0) = x implies thatx = δ᾽ A, sin(2n—1)xx/2l. Thus 
n=1 


2{' . (2n—1)xx 8i(—1)"*! 
A= 7 [| xsin θθ θα ἀπ SS τα 
8] So (--1) (2.---ἰ)πυε. (2n—1)nx 
and u= "2 ky (2n—1) cos y OF 


u satisfies the boundary conditions by construction, but as it is not 
permissible to differentiate twice under the summation sign though the 
separate terms satisfy the wave equation, it is not possible to show the 
series also satisfies the equation. u can be written as F(x —vt)+F(x+vt) 
where 


41 = (—1)"*! | (2n—1)nx 
F(x) = md, nee 3 


F(x) is related to the function sketched in Fig. 8 and is in fact the sine 
series of period 4/ representing x in [0, Π. By analogy with Fig. 8 it has 
a discontinuous derivative at x = 4I,4l,..., and its second derivative 
does not exist at these points. Thus, as any twice differentiable function 
of x+vt satisfies the wave equation, F(x+vt) are solutions of the wave 
equation for a given t at all but a finite number of points. 

The solution obtained is a solution in a somewhat generalized sense 
and in ihe general nomenclature of hyperbolic equations is called a 
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reflects the discontinuity inherent in the initial conditions. | oO 


= 


Simplified procedure In both Problems 2.1 and 2.3 the equation for 


X reduced to the simple harmonic one and it can be seen from these 


problems that this equation will always be obtained when the variable 
x only occurs explicitly in the equation as a constant multiple of 6?u/éx?. 
The same situation will also occur if there is a constant multiple of u in— 
the equation. Also, by analogy with the previous problems, possible 
solutions X can be obtained for those cases when u or du/éx vanish at 


« ῃ 
[ a 


εἰ 
ΓΝ 


‘weak’ solution. The discontinuity in the derivative of F samen 


we ., 


ν 


x = 0,c. The forms for two cases have already been found and the results 


for the other cases are 
(i) du/ex = 0,x =0,c: Χ = cosnnx/c,n = 0,1,2,..., 
(ἢ Gu/ex = 0.x τε = Ὁ x = δ; 
X = cos(2n—1)nx/2c, n = 1,2,.... 
Thus for equations of the form of (2.2) with a, = 0 and a,, a, constant 
it is possible to predict, for certain boundary conditions, without going 
through the general separation of variables procedure, the general form 


of the solution. It is useful to have the appropriate representations in the 


form of a table. 


Table 1. Representations for case a, = 0, a, ,a, constant 


Boundary conditions 


Hu=OQx=0,x=Cc 


(ii) du/€x = 0,x =0,x =c > ¥.(y)cos nnx/c 


(iii) u = O,x = 0; δι Οχ =0,x =c y Y (y)sin(2n — 1)xx/2c 
=1 


Ἂ 


(iv) du/éx = 0,x =O0;u=0,x =c Σ᾽ Y¥.(y)cos(2n— 1)nx/2c 


Entries (i) and (iii) respectively would have been appropriate for Problems 
2.1 and 2.3. Once the appropriate representation has been chosen from 
the table then substituting this in the relevant equation and equating the 
coefficients of the trigonometric functions to zero (cf. remarks in Problem 
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1.3) gives the equation for Y,. (Theorems 6 and 7 show that, provided u 
satisfies the relevant homogeneous conditions, the series can be dif- 
_ ferentiated term by term twice.) The approach is illustrated in the following 
problems. 


Problem 2.4 Solve the equation 67u/éx? = du/dt, O< x <1, t > 0, 
under the boundary conditions du(x,t)/éx = 0, x = 0; μά, ὃ = 0, t > 0; 
u(x, 0) = sin(zx/l),0 < x < |. (Physically this corresponds to calculating 
the temperature in the slab Ὁ < x < / when the end x = 0 is insulated 
whilst the end x = / is maintained at zero temperature, the initial tempera- 
ture being sin(zx/I).) 


Solution. Table 1(iv) gives the appropriate representation as u = 


EY T,(t) cos (2n — 1)x/2I. Substituting this into the equation gives 
n=1 
= (2n—1)? x? (2n—1)nx τὰ 41 (2n—1)nx 
προ Ὁ ρον ς-- 


n=1 n=1 
equating the coefficients of cos(2m— 1)zx/2l in this equation gives 
aT, (2n—1)* x? 
a 412 ly 
and the general solution is T, = A, exp[ —(2n—1)?x7t/4/7]. 

The boundary condition at t = 0 shows that Τ΄ (Ο) (= A,) are the 
Fourier coefficients in the expansion of sinzx/l as a cosine series of 
period 4/. This was considered in Problem 1.15 and A, can be found from 
that problem by replacing c by /. Hence 

ac yey | 2 
ΣῊ i: ν᾿ exp[ —(2n—1)°n*t/4l Ἰ ως 2 ax 
me (2n+ 1)(2n—3) 2l 
The coefficients in this series are exponentially damped for large n. 
Thus the order of differentiation and summation can be interchanged and 
u satisfies the differential equation; u also satisfies the boundary conditions 
(by construction). oO 


Problem 2.5 Obtain a solution u of Laplace’s equation finite everywhere 
in the sector 0 < r < a, 0 < 6 < a, where r,@ are polar coordinates, 
and such that it has zero normal derivative on the lines θ = 0 and 6 = ὰ 
and is equal to θ onr = a. 


Solution. u satisfies r? 0?u/ér? +r δι ὃν + 07u/d0? = 0, which is of the 
form of equation (2.2) with x replaced by θ. Table 1(ii) gives the appropriate 


representation for u to be δ᾽ F,(r)cosnz@/a. Substituting this in the 
a=0 
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equation and equating the coefficients of cos nz@/a to zero gives 
r?F" +rF’ —(n?n*/a*)F, = 0. 

For n # 0, independent solutions of this equation are r*"*’* and the 
finiteness condition shows that the positive sign must be chosen. For 
n = 0 the general solution is Alnr+B and for finiteness A = 0. The 
appropriate general solution for F, is thus A, r"** (all n) and the boundary 
condition on r = a shows that A, αὐτὰ are the Fourier coefficients of the 
cosine series of θ in [0, «], 1.6. 

1 {* 1 2.1" θ 
Ay = =| θ4θ---α, A,a™* = =| 8 cos "π΄ 4θ. 
α 2 α α 


0 0 


Integration shows that A, = 0 (n even) and A, , = —4a/n?(2r—1)?, 
[fr = 1,2,.. .] (cf. Problem 1.10(a)). The infinite series for the solution 
can now be written down. is 


The above analysis may be generalized to the case where u satisfies 
the equation (V?+k?)u = 0 where k is a real constant, and the boundary 
conditions are unchanged. The equation for Κ᾽, is 


αὐγὴ" +07 rF’ +(k?a7r? —n?n*)F, = 0, 


ΔῈ idF, n* 2 
» Ste get (1-S5)F. zi 


where z = kr. This is Bessel’s equation of order nz/a, and the solution 
finite for z = 0 is J,,,(z) which is the Bessel function of order ππία. Thus 


the solution of this second problem is found by replacing r”*/*, αὐτὰ in 
the solution found earlier by J,,, (kr), ᾽,,α α(Καὶ. 


nn/a 
Problem 2.6 Find a function satisfying Laplace’s equation, finite and 
continuous everywhere in the circle 0 < r < a, and equal to cos*@ on 
the circumference, r,@ being polar coordinates. 


Solution. There are no boundary conditions on θ = constant which 
would let us employ Table 1 but assuming a solution F(r)G(0) gives 


a 2 1d/( dF\ _ 
σ' = --λ6, 02) = AF 
G has the general solution A cos 1#@+Bsin j*@ and 4 can be found by 
noticing that G will not be continuous and single valued in [0,27] unless 


λ = πὖ (n any integer). The solution F, corresponding to 2 = n? can be 
obtained from Problem 2.5 on setting « = 7; the finiteness condition 
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vie 
Thows that F,, will be proportional to γῆ". Thus the general form for the 
Ἢ solution will be 


Ay ἘΣ γ"(Α, cos n@+B sin ηθ). 


Setting r = ashows that A_, B, can be calculated in terms of the Fourier 
coefficients of cos*@ and the latter are most easily calculated from cos 30 = 
4cos*@—3cos0. Hence B, = A, = 0, π #1,3,A,a=3, Α.α5 =1+ ΠῚ 


Problem 2.7 Find u satisfying V7u = δι δι, t > 0, finite everywhere 
_ within the sphere r < a, vanishing on r = a and equal to unity when 
t= 0. 

(This problem is equivalent to the determination of the temperature 
in a spherical region, the boundary of which is maintained at zero 
temperature, the initial temperature being prescribed.) 


Solution. The boundary conditions are spherically symmetric and it 
is reasonable to assume that u is a function of r only, hence 


1 ὃ 2 0u _ Ou 
Par\” or) at’ 
Writing ru as v gives 6*v/ér? = dv/ét, and in order that u be finite when 
+ = Oit is necessary for v to vanish at r = 0. v will also vanish on r = a 
2 and, as v satisfies an equation where r only occurs in the term 67v/dr . 
__ Table 1(i) shows that the appropriate representation for v is 
οο 
>> Τι(θϑίη ππνία, and substituting this in the equation gives 


n=1 

Τ' = (n’n*/a’)T,, with a general solution T, = A, exp(—n?n*t/a’). The 
condition at t = 0 shows that A, are the Fourier coefficients of the sine 
series of r on [0,a], ie. A, = (2/a) [* rsin(nnr/a) dr = 2(—1)"**a/nn. O 


Problem 2.8 Solve the equation 07u/éx?+6?u/dy? = é7u/dt?, t 0, 
0<x<c, 0< y<d, under the conditions u=0, x =0, x =c; 
y =0,y = d:u = xy(c—x)(d—y), δι δὶ = 0,t =0,0<x<c0<y<d. 
(A physical interpretation of this is the determination of the shape of 
a rectangular membrane given its initial shape and that its initial velocity 
is zero.) 
Solution. As this problem involves three independent variables Table 1 
is nee directly applicable, but assuming a product solution gives 
Σ᾿ Τ᾿ » dats al ae 
Ὥς ΤΑΣ τῷ ἂν 
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Both sides of this equation must be constant ie. 
xe we y eg 
—+—=-4 —=-i. 
ΤΠ 
As X and Y are functions of independent variables x and y, both the 
terms in the first of the latter equations must be separately constant, i.e. 
γ" eer x" 
Fa τ 
The equation for Y is similar to the equation of Problem 2.1 and in 
order that u vanishes at y = 0,d it is necessary that μ = m?x?/d? (man 
integer). Similarly from the equation for X and the conditionsatx =0.c 
it is found that (A— μὴ) = n?x7/c? (n an integer). Substituting A, μ in the 
equation for T shows that the general form for u is 


> y [A,,,cosK, t+B_ sin Καὶ ‘]sin—sin—, 


n=lm=1 a. 
where ΚΖ, = n?x?/c?+m?n?/d*. du/ét will vanish fort =O if B,, =0 
and the other condition at t = 0 gives ; 


= μ--λ. 


xy(e—x)(d—y) = δ he a sin— sin =, 0<x<c0<yK<d. 
This is a double Fourier series and, for a general left-hand side, A, 
could be found by successive application of ve (1.5), ie. 


mo aI I. xy(e—x)(d— — y)sin— ~ sin = dxdy. 


This elaborate sree can be avoided in on case as the integrand is a 
product of a function of x and a function of y. The sine series for x(c—x) 
has been obtained in Problem 2.1 and that for y(d—y) can be derived 
from this by replacing x,c by y,d. Hence A, = 0 if either m and/or πὸ 


is even and 
64c7d? 
Oa = n°(2r — 1)>(2s—1)>” rs= | ae Se 
Hence 
«ἜΣ Ὁ ὦ cos Boa COS K>,_ 125-1! Oe ann 1590-4: is (2r— (2r—I)ny ς in US yx O 
δ 2 (2r — 1)5(25 -- 1)" pail al ς ; 
Problem 2.9 Solve the equation 
Ctu «1 070 
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where a is constant, under the conditions u = 67u/éx? = 0,x τ 0,x = |, 
—t > 0; u(x, 0) = f(x), du(x,t) δὲ = 0,t = 0,0 <x <1. 

(This can be interpreted physically as determining the transverse 
displacement of an elastic beam of length /, simply supported at its ends. 
The initial velocity of the beam is zero and its displacement is f(x).) 


Solution. The equation is not of the form of (2.2) but assuming a product 
solution X(x)T(t) gives X"”’—AX = 0, T''+Aa*T = 0, where / is the 
separation constant. Making the trial substitution X = Ae” gives 
μ᾽ = d,ie. μ = +/+, +id* hence the general solution is 

X = Acosj*x+B sin A*x +C cosh 4*x + D sinh A*x. 
The conditions at x = 0 give A+C = 0 = A—C; ie. A = C = 0; the 
conditions at x = ἰ give 


Bsin λὶ! -- Γ sinh 2*] = 0 = Dsinh 4*!—Bsin #1. 


Thus, if B and D are not to vanish, either sinh A*/ = 0 or sin A*/ = 0. 
The former gives ἃ = Ὁ (purely imaginary values for j* are roots of 
sin 4*] = 0) which implies X = 0; hence 4* = nn/l where n is any 
integer. Hence D = 0, and the appropriate representation for uis thus 


Σ᾽ 7,(t)sin nxx/l where Τῇ +(n*x*a*/I*) T, = 0.The general solution for 
n=1 

T, is A,cosn?n?a*t/I?+B,sinn?n?a*t/I?, δι δὲ =0 at t =0 implies 
B, = 0 and u(0,t) = f(x) implies that A, are the Fourier coefficients of 
the sine series for f(x) in [0, /], 1.6. 


I 
ἡ μὰς | f(x)sin—— dx, 
fe l 
Thus, for any given f(x), a complete series solution has been obtained. 0 


2.3 Inhomogeneous Problems It is necessary, in order to solve equation 
(2.1) for f # 0 with u or Gu/dx taking prescribed but non-zero values on 
x = 0,c, to modify the techniques used in the previous problems. 

In some cases f and/or the form of the boundary conditions on x = 0,c 
are sufficiently simple for a function U (often a polynomial in x) to be 
found by inspection, which satisfies both the equations and the conditions 
on x = 0,c. It is unlikely that a simple function U can be found satisfying 
whatever conditions are prescribed on y = constant also, but writing the 
unknown u as U+p will give an equation of the form (2.2) for v with v 
satisfying homogeneous conditions on x = 0, ες. The problem for v will 
be of the general type discussed earlier. 
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Homogeneous boundary conditions For problems where u or du/éx are 
prescribed to vanish for x = 0,c (but it is not possible to use the above 
simple approach), it is only necessary to modify the previous approach 
very slightly. Ifthe problem forf = Ocan be solved by usinga representation 
for u as a Fourier series, then a useful method to try for the inhomogeneous 
problem is to assume a Fourier series in x of the same type for f Sub- 
stituting both forms into the relevant partial differential equation and 
equating coefficients of trigonometric functions enables a differential 
equation to be obtained for the Fourier coefficients. Both the above 
methods are illustrated in the following problems. 

Problem 2.10 Find a function u satisfying the equation and conditions 
of Problem 2.3 except that at x = / which is replaced by @u/dx = 1. 


Solution. x satisfies both the differential equation and the conditions at 
x = 0,1]. Thus writing u as x +v shows that v satisfies the same differential 
equation, and the boundary conditions on v at x = 0 and x = / are those 
imposed in Problem 2.3. The conditions u = du/dt = 0 at t = 0 give 


v = —x,év/6t = 0 at t = 0. Hence —v is the function obtained in 
Problem 2.3 and 
δ᾽ τὸ" {-- 752 πυι. πΧ 
“τε ὩΣ eee gpa ap oO 


Problem 2.11 Solve 07u/éx?+é07u/dy? = 2, O<x<c, O< y<d, 
under the boundary conditions u = 0 on x = 0,c:0 « y<d:u=0O, 
y=Qcu/oy=Oy=¢0<x< 


Solution. The right hand side of the equation is sufficiently simple to 
suggest seeking elementary solutions. Particular solutions are x* and 
γ, and adding to either linear combinations of x and y will still produce 
solutions of the differential equation. If we consider x? +ax+b then the 
conditions at x = 0 and x = ς are satisfied by b = 0, a = —c. Hence 
writing u as x*—cx-+v shows that v is a solution of Laplace’s equation 
vanishing on x = 0,c;the conditions onu at y = 0,dshow that v = cx—x’, 
y = 0; δυῶν = 0, y = d. The appropriate representation for v is that of 
Problem 2.1, i.e. 


v= 4 [A, cosh nzy/c +B, sinh nzy/c]sin nnx/c 
n=1 
The conditions on y = d will be satisfied if A, sinh nnd/c = — B,coshnnd/c 
and that on y = 0 shows that A) are the coefficients of the sine series of 
x(c—x) in [0,c] and thus A, have the same values as in Problem 2.1(ii). 
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δὰ a+ te 
An alternative approach would be to satisfy the boundary conditions 

on y =0 and y = d; y?+ay+b satisfies both the equation and these 
_ conditions ifa = —2d and b = 0. Hence writing u as y?—2dy+w shows 
4 _ that w satisfies Laplace’s equation and the appropriate homogeneous 
_ conditions on y = 0,d. On x = 0 and x = οὖν is equal to 2dy—y’. 
_ The boundary value problem for w can be solved by means of a representa- 

tion of the form of Table 1 (iii). oO 


Problem 2.12 Solve 0?u/éx*—6?u/dt? = exp(—1), t > 0, 0<x <1. 


τς under the conditions μ(0, 1) = u(l,t) = 0, t > 0; u(x,t) = Au(x, t)/et = 0, 
 t=0,0<x<l. 


4 a? Solution. It is not easy to find by inspection a function satisfying both 
__ the equation and the boundary conditions and we need to use the more 
| ΗῚ general approach described. If the right-hand side of the equation were 


: zero then the appropriate representation would be u = x T(t) sin nzx/I 


᾿ς and we therefore need to obtain a sine series representation for exp( s t) in 
[0, Π. This representation will of course be exp(—1) times the sine series for 
unity, the latter series can be calculated by the methods of ὃ 1.2, and it is 
found that 
aa 


b 


θ τ χ ce. 


WY 4 y oe eee 
m ou 2r—1 


Thus substituting this series and that for vin the above differential equation 


gives 


a - ΠΥ ἀφ ΒΡ ΝΣ ππχ 4 © sin(2r—1)nx/l 
= ——,- T. -—T” | sin— = - --------- -- 
ἫΝ >| Pon *n | ΠῚ ἊΣ Bat τ τς, αὐ 
: Equating coefficients of the trigonometric functions gives 
a. 4γ32π2 
": — a IG, =O Pm Ὁ.1.2 τ, 
~  (2r—1)?n? 4εχρί-- ἢ 
"a ——_,—_T,_,-T_, = —"—- r= 1,2, 
δὰ Ι2 2r-1 2r-1 r—1)’ RSENS ie 
ἊΜ m(2r—1) 
"ἢ 


ΝΜ ᾿ 7 
. ; Ἢ Ἦ 4 » Ἢ 
beat ἊΝ ἡ ‘eth ΥΎΝ ΡΝ ΤΥ le ν᾿ a” : 
i] - γ᾽ ay 4 . es be, y » ὙΠ. Ἅ. λ } : ΕἾ = , Ἴ 
ΨΥ aa, ΡΥ | = , _ * g ; ᾿ J i 
“s » δε a 4 [= γ .» ι ὦ a" 


Πα = du/ét = Oat t = O imply that T, = Τ' 


follows that T,, = 0. The general solution for T,__, is 
4]? exp(—t) (2γ --- 1)πὶ 


νοι = παρ  ΉΓγ-- ΠΕΣ +P] 


+A, _, COS 
+B,,_, sin 


The conditions at t = 0 give 
A ἐλ 412 
ὅτι πίΖε-- 1) [ 1? +(2r — 1)2π27᾿ 


4|" 


Ῥ ~ 222 1)?[1? + (2r—1)?x?] 


2r-1 
Substitution of these values into the expression for T,,_, gives 


412 <= 1 (2r—1)xt 
sity ae (2r—D [FP +(2r—1)*x"] expt —1)-c0s 


=) | i 


Problem 2.13 Solve 07u/éx* —du/ét = exp(—t)sinx, 0 < x <c¢,t >0, 
under the conditions u(0, t) = 0, du(x,t)/éx = 0, x = c;t > 0, 
u(x,0) = 0,0 <x <e. 


n=1 
for sin x. This is obtained as in § 1.3, and 
ane. Ἢ (—1)"*! __ (2n—1)nx 
sinx = 8ccosc bs [(Qn—1pn2—4c7] 5ς ᾿ 
n=1 


in this case the representation of sin x would be itself. ] 
Substituting both forms in the equation gives 


κι [  (2n-1)?n?_ 47. -(2n—1)nx | 
ΤᾺ 2: ᾿ς Te ay = 8ccoscexp(—t) 


τι. (-1* Geet 
pa Qn—1Pn—4c? ΟΣ 
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ΠΑΡ ΡΥ 

5 . . ads Ψῃ τὴ 
= 0, for 81} η when t = 0, and ᾽ 
as the equation and the conditions for T,, are both homogeneous ἡ 


Σ᾽ Τί sin(2n— 1)xx/2c, and we therefore require a similar representation 


Ὶ 
γὴν 


(2γ -- 1)πὶ a 
an 


l . (2r— ἐπε sin ΒΈΤΠΕΣ, oO 


ι 
ig \ ey 
be i 


Solution. If the right hand side of the above equation were zero itfollows 
from Table 1(iii) that an appropriate representation for u would be 


[It has been assumed that c is not equal to (2n—1)n/2 for any integer n; Ὁ 


a4 ᾿ 


and hence 
aT, . (2n—1)?x? _ 8ecosc(—1)"exp(—t) 
dt 4c? κ΄ (2n—1)?n?—4c?. " 
The general solution of this equation is 
32c? cos c(—1)"exp(—t) 
[ (2n—1)?n* -- 4.31: 
and to satisfy the condition on t = 0 we require that 
debe A 32c* cos εἰ -- 1}" 
n [(2n wh. 1)2π2 -- 4.272 μ 


T, = A,exp[ —(2n—1)?27t/4c?]+ 


Hence a complete series solution for u may be written down. oO 


Problem 2.14 Solve 07u/0x? = 6du/ét,0 < x < |, t > 0, under the con- 
ditions Ou(x, t)/6x = 4, where q is constant, for x = 0, u(l,t) = T, where 
T, is another constant, and u(x,0) = 0. 


Solution. Any linear function of x satisfies the equation and the particular 
choice q(x — ἢ Ὁ T, also satisfies the boundary conditions. Writing u as 
q(x—l)+ T,+v shows that v satisfies the same equation as u, and also 
satisfies homogeneous condition for x = 0, and x = +/ and 
u(x,0) = -- Τὶ -- αίχ -- ἢ. The appropriate representation to choose for v 
will be that of Problem 2.4 and the results can be obtained from those of 
that problem by replacing sin nx/l by — Τὸ — q(x -- ἢ); this means that 

v= Σ Α, οχρί —(2n— 1)?x1/4P Joos ΠΕ 
n=1 


where 


4(— 4(—1)"T, 8lq 
(Qn—1)n -- Ι)ὴπ (2n—1)'n? ἢ 


Inhomogeneous boundary conditions When the prescribed values Ἢ 
u or Gu/Ox on x = 0,c are not zero and the problem cannot be reduced 
to homogeneous form, then a general, though standard, approach has 
to be used. If the corresponding homogeneous problem could be solved 
by using a particular Fourier series representation then the technique is to 
multiply the differential equation by the general member of the Fourier 
series [e.g. sin nmx/c for the case of Table 1(i)] and integrate the equation 
with respect to x from x = 0 to x = c. Integration by parts of the term 
involving 67u/éx*, [for the case when Fourier series are applicable there 
will be no term involving du/éx], will lead to a differential equation for the 
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l 
A,= -7| [T, Ὁ (χ -- ἢ] οοβ(2η -- 1)πχ 2] dx = 
0 


ΠΎΙ, ΜΕ Ae SP Υ δ ee” Pee Oe eye |! ΨΚ Ὰ 


-: 


4 


Fourier coefficients of u. This technique is illustrated in the following 
problems. 


Problem 2.15 Solve the equation of the previous problem under the 
conditions 
u(O, t) = t, u(c,t) =0,t > 0, ux,0) =0,0<x<c. 
Solution. The corresponding homogeneous problem [i.e. 
μ(0, 1) = u(c,t) = 0] could be solved by a sine series representation and 


we thus multiply both sides of the differential equation by sin nzx/c and 
integrate with respect to x from 0 to ¢, 1.6. 


2 
| Ὁ τ βίη 5 8 ΡῈ ἐὰ -| ἢ oe sin ἀκ. 


Integration by parts twice of es left-hand side of this gst: gives 


Gu... MEX Pr". RE nx =" n?x? nt 
— sin —— -- ---ἰ ucos—— ie ὩΣ 
ΟΧ Pau, © se roan c 


x=0 


@ 
u(x, t) can be represented in0 < x < cas δ᾽ b/(t)sinnnx/c, where 


᾿ : n=1 
b ΞΞ - | u(x, t)sin—— dx. 
oA FS c 
Thus, on applying the conditions at x = 0,c to the equation obtained 
by integration of parts, 
db, nn? πὶ 
@ eu | aw 
i.e. the required equation for b, . The general solution is 


5 ΨΕΥ 2 
b, = A,exp(—n*n*t/c )+2( + “| 
and in order to satisfy u(x,0) = 0, A, = 2c?/n?x° giving 


i= ) ee [exp(—n 22t/e2)— 1]sin™ + = 2, “sin. 
n=1 a= 
The last series is very slowly convergent and for computational pur- 
poses is not satisfactory, but its exact form is given in equation (1) Appen- 
dix 1 and use of this result and (7) of the same appendix considerably eases 
numerical computation. In general the results of Appendix 1 are often 
useful in avoiding numerical calculation of slowly convergent series. 
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It would not have been correct in this case to have differentiated the 


sine series twice to calculate the Fourier series of 67u/@x? as the sine series 
only represents u in the open interval and cannot be differentiated term 
by term [c.f. remarks preceding Theorem 6]. An alternative method would 
have been to use Theorem 6(ii) to calculate the cosine series of éu/dx and 
then applied Theorem 6(ii) to this series to obtain the Fourier series for 
67u/dx*. The above technique involving integration by parts is more useful 
in that it avoids the necessity for remembering formulae. 

It would also have been possible, with some ingenuity, to have reduced 
the problem for u to a homogeneous one. A simple function satisfying the 


boundary condition is u, = t(1—x/c), but 


a ae 
and it is therefore necessary to add to u, a function f(x) such that 
f” = —1+x/c and with f(0) = 0 = f(c); such a function is 
(x? —3cx? + 2c?x)/6c. Hence writing u as 


l 
(1 ~*) +260 —3er +2c*x)+v 


yields a homogeneous problem for v. The form in which u is obtained in 
this case is precisely that obtained by using the results of Appendix 1 to 
transform the series solution obtained previously. oO 


Problem 2.16 Solve 0*u/éx* = 67u/ét?, t > 0, 0 < x < c, subject to 
the conditions u(0,t) = 0, du(x, t)/dx = sinwt,x = c,t > 0: 

u(x,0) = du(x,t)/et = 0, t= 0, O< x <c [ὦ is not an odd integer 
multiple of /2c]. 


Solution. The appropriate representation for the homogeneous problem 
would be that of Table I (iii) and we thus multiply both sides of the equation 
by sin (2n—1)zx/2c and integrate with respect to x from x = 0 tox = ¢, 
Le. 


(2n— jax | (2n—1)nx 
ἃ: πο ax = | oe 0 a 
Integration by parts twice gives 
Ou __ (2n—1)ax _ (2n-1)n «τος 21- πα 
ΟΧ ae. 7 ΕΑ 2c 2c Mes 
(2n—1)?n? (© | (2n—1)xx a [Ξ . (2n—1)xx 
4c? : usin—~- dx = af J) a, Pa dx 
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If b, is defined ‘by u= ? b (t)sin(2n—1)xx/2c it follows, on using the 
prescribed conditions, that 

d*b. , Qn- 1)?x*b, a A= 1p" sin wt 

= ag 4c? ς : 


The solution of this satisfying the conditions that db,/dt = b, = 0 at 
t = Ois 


84 — 17"! E Jah 2ως in on 


(2n—1)'2? —4c*w? (2n— i 2c 
Thus the series solution for u is obtained. oO 
Problem 2.17 Solve the equation 


ay Ou δὲ 
dx* αν dy* 


under the conditions u(0, y) = u(c, y) = 0,(07u/dx?), _. = (07u/dx’),_. = 
2,0 < y < δ᾽ u{x,0) = u(x,d) = 0, du/oy =~ 0, y = 0,42,0<x<c. 
(This problem can be interpreted physically as calculating the deflection 
of an elastic plate, two of whose edges are clamped and the other two are 
simply supported. It is also acted upon by bending moments uniformly 


distributed along the free edges.) 


=0 0<x<c0<y<d, 


Solution. This boundary value problem is of a type not previously 


considered and it can be verified that the method of separation of variables 
is not suitable. The problem can be somewhat simplified by making all 
the boundary conditions on x = 0 and x = c homogeneous and this can 
be achieved by setting u = x(x—c)+v and v will satisfy the same equation 
as u. The conditions on v show that it can be represented in [0,c] by a 


series of the form δ᾽ Y,(y)sinnmx/c and, as 07v/@x? vanishes at x = 0, ¢, 
n=1 
it can be shown from Theorem 6 that it can be represented as 


— ¥ (n’n*/c’)Y, sin nzx/c for x in [0,c]. Further application of Theorem 


n=1 
6 shows that the series representing 0*v/dx* is δ᾽ (n*x*/c*)Y, sinnzx/c. 
n=1 
Substituting the series in the equation and equating coefficients of 
sin nzx/c gives 
4n*Y —2n?n*c?7¥"'+c*¥,"" = 0, 

the conditions on Y, are that Y’ = 0 at y = 0,d and Y,(0) and Υ (4) are 
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os f oe » y oe ay. 
΄ wt De 1 ἴ pe ΟΣ d 
᾿ τω» ῃ ες 1 ὺ ᾿ ὃ 1 ᾿ as st 


the ΣΡ αἰποίοηῖς of the Fourier sine series of χίς -- x). Tistc have been oe 
Ἴβηρα in Problem 2.1 and thus Υ(0) = Y(d) = 0 for n even and 
. 8c? 

. a | Y,,- 10) = Y,,,(@) = Ὡν-- ηἾπϑ᾽ = 1.2... 
“Τῆς general solution for Υ is a linear combination of cosh nmy/c, 
sinh nzy/c, ycosnzy/c, ysin nny/c and there will be four conditions to 
. ᾿ determine the four constants. The calculations can be simplified by notic- 
: _ ing that Y, satisfies the same condition at both boundaries, i.e. it is sym- 
; ᾿ metric about y = 4d and we choose the linear combination with this 

aha 1.6. 
Υ - A, cosh — Pe —4d)+B(y— 4d) sinh — (y— 4d). 

δ τὰ it follows from ts conditions ¥(0) = Y¥(d) = Y/(0) = Y/(d)=0 

that A, = Β -- 0, thus 

) (2r—1)nd δι (r~Dad > Be 

| A,,_, cosh ze +8,._, 5 sinh ς = Groin 
) ω- 1)x fae: inh (2 - Ln (2r —1)nd 
δ Ὁ 2c 2c 


, va = 
ἫΝ a og ate ἃ 


+B,,_, sinh 


| ALY 


_ These equations can now be solved to give A,,_,, B,,_, and the series 
solution for v is then found. 0 


EXERCISES 


md, Solve 67u/@x? + @7u/dy? = 0,0 < x < a,0 < y < b, under the con- 
ditions u(0, y) = μία, y) = 0, 0< y «5; u(x,0) = 0, du/éy = 4, y = hb, 
Ὁ <x <a, where q is a constant. 


2 Obtain a solution u(p,@) of Laplace’s equation finite everywhere 
within the circle of radius a, (p, 0) being the usual polar coordinates, for 
_ the two cases 


(i) u(a, 0) = sin @, 


(ii) μία, 0) τ, οἰὶηθ.0 «θ «π 
= 0,2 < θ᾽ ες 22. 


+ 3. Solve the equation of exercise 1 under the conditions u(O, y) = 0, 
— διμδχ = 0, x = a,0 < y < b; u(x,0) = 0, u(x, δ) = 1lO<x<a. 


_ 4. Solve the boundary value problem of exercise 1 when u is prescribed 
to be unity on x = 0,x =aand y = 0. 


= a ΕΣ Β 


Νὰ 


he bat , fide ἀν; : | ἡ δ Ὑ. es a? i“ : ἜΝ ies 

, εἰ . La ἥν "EF ἢν, 
ἐν ἢ Solve ae ὌΧ 0 «χε τῇ t>0, ‘subject to the c condita : 
u(0,t) = a, u(t) = b, t > 0; u(x,0) = 0, O< x <1; a and δ᾽ cing — 
arbitrary constants. fest A 


6. Solve 67u/ét? = 07u/dx?,0 < x < lI, t > 0, subject to the conditions — ; 
u(0,t) = 0, δι Οχ τὸ 1, x=], t>0; u(x,0)=0, dufot=0, t=0, © 
Oe εν a 
ie 
7. Solve the previous exercise when the condition on x = / is u = sint 
[nx # | for any integer n]. 


af 


᾿ 
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᾿ Fourier Integrals 


where 


Ψ͵ 
"ὃν 


3.1 The General Fourier Integral Ifa function f(x) satisfies the conditions 

of Theorem | for all intervals then (by setting a = —c in the expressions 
of 1.1) a Fourier series representation for f(x) can be obtained in [ —c,c] 
for all values of c. This series will not, however large the value chosen for 


τς ὦ, represent f(x) for all x unless it has period 2c. We therefore examine 


the possibility of extending the ideas of ὃ 1.1 to represent a function for 
all x in terms of trigonometric functions. 

In the following analysis it is convenient to use the exponential form 
of the Fourier series described immediately prior to the statement of 


4 _ Theorem 1. Any function f(x) satisfying the conditions of Theorem 1 is 


such that, except at a discontinuity, 


ἀν * 


a= - ρος 
c 


1 
“= =| f (x)exp(—innx/c) dx 


a,exp(innx/c), |x| <c, 


where ch = 7. 


-- i | ; f (x)exp(—inhx) dx, 
> a (a 


_ Hence f(x) = = : hF(nh)exp(inhx) 
FB) = |" f(x)exp(—ifx) dx. 


The ordinary definition of an integral as a limit of a sum now suggests 
that, on taking the limit ἢ > ty 


70) =— =| exp(iBx)F(B) dB (3.1) 
F(B) = lim F(f) = ᾿ f (x)exp(—iBx) dx. (3.2) 


F(B) is known as the Fourier transform of f. 
The above analysis is clearly not rigorous and assumes existence of 


where 


- various limits and integrals; equations (3.1) and (3.2) are, however, valid 


under fairly general conditions, and their validity can be formalized in 
the following theorem. 

Theorem 8 If f(x) is piecewise smooth on every finite interval of the 
x-axis and {* | f(x)| dx exists, then the right hand side of equation (3.1) 
_ represents 
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¢ 
f 


(i) f(x) at all points at which it is continuous, 

(ii) ΣΓ f(x, +0)+ f(x,—0)] at any point of discontinuity x = x,. 

An important result "associated with Fourier integrals is the Riemann—_ 
Lebesgue lemma which states that, if f satisfies the conditions of Theorem 
8, lim F(f) = 0. 

Bom 

The Fourier integral is of such general applicability that tables of 
functions and their corresponding transforms have been compiled; 
some relevant results are tabulated in Appendix 2. Thus only very few — 
examples of direct calculation of Fourier transforms will be given. The 
process of determining f(x) given F(f) is known as the inversion of a 
transform and this process is sometimes helped by the following theorem. 


Theorem 9 (Convolution Theorem) If f and g are functions satisfying the 
conditions of Theorem 8 and F(f) and G(f) denote their Fourier trans- 
forms then the inverse of FG is{”  f(t)g(x—t) dt. 


Problem 3.1 Find the Fourier transform of the function f(x): f(x) = 
x <0; f(x) = exp(—x), x > 0; f(0) = 4. Hence obtain an integral 
representation of f(x) for all x. 


Solution. From equation (3.2) it follows that, since f...0 for x < 0, 


oak ais ge pola ie 
ro =| exp(—x—iBx) dx = [+ip ~ 1467 


f(x) satisfies the conditions of Theorem 8, and hence 


f(x) = exp(ifx)F(B) dB = -- : a: τ τ exp(iPx) dB 
* (1—iB)exp(iBx) dp n rm — ip)exp(iBx) 4 
= 5, 1+ p? _» I+f (ἃ: 
B can be base's by -β τ the second integral thus: 
Fic ” [cos Bx +B sin Bx] 
1+ p? is - 


Problem 3.2 Find the Fourier transform of exp(—a|x|) (a Ὁ 0) and 
hence obtain an integral representation for exp(—a|x|). 


Solution. 


F(B) Ἢ exp(—a|x|—iBx) dx 


δ exp(—ax—iBx) dx+ [ ᾿ς exp(ax -- ἰβχ) dx 
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a+ip α--ἰβ a*+f? 
ΝΣ feces Theorem 5: 


exp(—a|x|) = a0 See [] 


The above identity is not very easy to obtain by elementary means, but 
one method would be by means of contour integration and use of Cauchy’s 
theorem and readers familiar with complex variable theory might be 
interested in this alternative approach. [exp(ixz)/a?+z?] is integrated 
over the closed contour in the z plane consisting of the real axis from 
z= —R to z = R and the semicircle in the upper half plane joining 
these points. There is a simple pole at z = ai with residue —iexp(—ax)/2a, 
and Jordan’s lemma shows that the integral over the semicircle tends to 
zero as Καὶ — οὐ, thus applying Cauchy’s theorem and taking the limit as 
R — οὐ gives the result for x > 0. The corresponding calculation for 
x < 0 requires the contour to be in the lower half plane; the case x = 0 
15 elementary. 


_ Problem 3.3 If f(x) and f'(x) are continuous functions of x satisfying 
all the conditions of Theorem 8, show that the Fourier transform of f' 
is 1β F(B), where F is the Fourier transform of f. 


Solution. The required Fourier transform is 
7" Πρλεχρί- ἰβχ) dx = [f(xexp(—ipx)] τ ὃ 


+ip rr f(x)exp(—iBx) dx, [on integration by parts]. 


Since the integral of | f(x)| from x = —c to x = οὐ converges it is 
necessary that f 0 as x + +00; thus the integrated terms vanish and 
the required result is proved. 0 


Problem 3.4 Show that an even function f(x) satisfying the condition 
of Theorem 8 has a Fourier representation of the form je G(B) cos Bx dp. 


Solution. The Fourier transform of f(x) is 
F(B) = |", f(x)exp(—iBx) dx = 
οὶ [: “ £(x)exp(—iBx) dx+ ἫΝ °  f(x)exp(—iBx) dx. 
x can be replaced in the second integral by —x and the evenness property 


then gives F(f) = 2 J5. f(x) cos Bx dx. Cos Bx is an even function of β 
hence the same is true of F; a particular example of this general result has 
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been derived in Problem 3.2. f(x) is such that Theorem 8 can be applied 
and thus 


fix) = =|" Fexptips) dp 


= = % F(B)exp(iBx) dB+ - ee F(B)exp(iBx) dB. 


Replacing β by — f in the second integral and using the evenness property 

of F(B) gives f(x) = (1,πὴ [τ΄ F(B)cos Bx dB, which is of the required 

form. oO 
It can similarly be shown that, for f an odd function of x, 


f(x) = (1,πὴ ia G(B)sin Bx dB where G(f) = 2 = f (x)sin Bx dx. 


3.2 Fourier Cosine and Sine Transforms By constructing even and odd 


extensions of a function f(x) defined for x > 0, one can (from Problem 
3.4) obtain integral representations involving only sines or cosines | this 
is analogous to obtaining sine and cosine series representations on a 
finite interval]. Hence we have the following theorem. 


Theorem 10 For a function f(x) satisfying the conditions of Theorem 8 
for x > 0 the Fourier cosine and sine transforms F (f) and F (f) defined by 


FB) = |” f(x)oos Bx dx, (3.3) 
F(B) = |” f(x)sin Bx dx (3.4) 
are such that at all points x > 0 at which f is continuous, 
f(x) = 2/n) | F.(B)cos Bx dp, (3.5) 
f(x) = (2/n) |" FB)sin Bx dB. (3.6) 


At those points x = x, at which f is discontinuous the left hand sides of 
equation (3.5) and (3.6) have to be replaced by 2 f(x,+0)+ /(x,—0)]. At 
x = 0 the right-hand side of equation (3.5) converges to {(0+0) and that 
of equation (3.6) to zero. 


oo 


Problem 3.5 Calculate the Fourier cosine and sine transforms of a 


exp(—x)cos x. 


Solution. 
F (B) = [ξ exp(—x)cos x cos Bx dx, F,(B) = |* exp(—x)cos~x sin Bx dx. 
The most direct method of calculating these integrals is by expressing 


45 


aed 


i * Ρ 


ot 
AS rans ahd sin Bx in terms of exp(+ifx) and taking cos x as the real part 
Ε οἵ ae) Thus, on defining G(B) by 


by: 


πα... 
ΝΎ: 
ν᾽ 
[ὦ 


F(B) = 


exp(—x)sinx = Ξ 


|, f(B)cos Bx dB = 


2 Re[G(B) + G(—)], 


G(p) = [exp —x+ix+iBx) dx = 


we have that 


Ε(β) = 


Hence 


F(p) = 4Re —i[G(6)—G(—f)]. 


p*+4 


Ὁ (B* +2) 
p*+4 


1 
1—i(B +1) + 


_ 1Re2(1+i)(2+ p?)—4ip? εἰ 


(2+ B?)? --4β2 
1 Re 4i(1+i)B+2f(2+B?) _ 


© B® sin Bx 


_ 146 +1) 
1) ~ 242p +B?’ 


B? +2 
pt+4’ 
B? 


᾿ς β-4 


From these expressions and Theorem 10 are obtained the identities 


exp(—x)cosx = 4 i, cos Bx dB 


dp. 
p*+4 
Problem 3.6 Solve the integral ΡΝ 


Solution. Theorem 10 effectively states that if 


then 


Problem 3.7 Find the Fourier cosine transform of exp(—x?). 


J> cos Bx exp(—x?) dx. The simplest method 
of evaluating this integral is by a complex variable method but we shall 
_ consider an alternative method which uses more elementary ideas. The 


Solution. We have F (8) = 


\° f(B)cos Bx dB = F(x) 


ve: x < I, 
ἢ x > 2b 


| ΠῚ = (2/n) 7" F(x)cos Bx dx, 
and applying this result gives 


τῳ 


== | cospxdx = 
Ὁ sin B 

δος, ἃ 

whine 


2 sin B 


πβ 


β Substituting this into the original equation and setting x = 0 shows that 


O 


0 and thus ΔΕΒ can be evalu. τῇ 


integral is uniformly convergent for x > 

ated by differentiating under the integral sign, giving it =" 

dF./dB = — Ε x sin Bx εχρί-- ἢ) dx = 4 \5 sin Bxd exp(— x?) a 
-Ὁ}Ρ; Εν. 
2 4 


(on integration ὃν parts). a 

Hence F, = A exp(— —1?) where A is a constant. A may be found by “-. 
setting β = 0 and using the known result 2 [Ὁ exp(—x*) dx = Jn. The 
final saute is 2Ρ, (β) = «πεχρί-- 23). ΓΙ ΣΝ 


[For readers familiar with complex variable theory the following 
direct method might be of interest. As e~*’ is an even function of x it 
follows, by an analysis exactly similar to that of Problem 3.4, that 


2F (B) -}" exp(—iBx—x*) dx = οχρί-- ee exp(—x+}ip)? dx. i 


exp(—z’) is an analytic function of z and thus, integrating it round the a 
rectangular region with vertices at z = +R, +R+4if, gives ' 


sip x 
at exp(—z?) dz+ ἶ ἫΝ" exp(—z”) dz+ Ἢ ag exp(—z?) dz+ Γ ᾿ 
- τ εἶλΡίτ- 2) dz = 


As R — co the second and fourth integrals tend to zero, and setting — 


= x+4if in the third integral gives iS 

Ἦ exp(—z’) dz = τὰν exp[ —(x+4if)*] dx. ag 
The left-hand side of this identity is known to be «(π and thus F(8) can _ 
be calculated. | oO ἊΝ 


Problem 3.8 If f, f’ and f” are continuous functions of x satisfying all 6 
the conditions of Theorem 10 determine the cosine and sine transforms 22 
of f” in terms of the corresponding transforms of f and the values of f | 
and f’ atx = 0. ως 
Solution. The cosine transform of f” is Ο,(β) defined by ne 
 G,(B) = \ f’cos Bx dx ς a 
and integrating by parts we get 
G(B) = [f'cos Bx} τ +B : f'sin Bx dx 

= [f'cos Bx]""” + BLS sin Bx te a f cos Bx dx 

(on a further integration by parts). f, ζ΄ both vanish as x > oo and thus 
ay) 


} G.(B) = —[/f'],..—B’F (8). Repeated integration by parts of the cor- 
responding integral for G.(f), the sine transform of f, gives 


G(B) = BLS], - o ~ B’F (8). ΕἸ 
3.3 Solution of Boundary-Value Problems Using Fourier Transforms 


Fourier cosine and sine transforms are useful for the solution of boundary 
value problems for 


07u 07u du 
Bega TOs +h b)s + hsb = fayx>0<y<d, 


_ where a,, a, are constant, with u - 0 as x > οὐ and either u or du/éx 


_ prescribed on x = 0. The general method is to take either the sine or 
cosine transform of the above equation; this gives (using the results of 
Problem 3.8) 


a,[—B?U,—(éu/éx), _ 5]+a;U,+ [μι cos Bx dx = F., (3.7) 
a,[—f°U,+fu),_ o]+a,U,+ |"LusinBxdx=F, (38) 


where capitals denote Fourier transforms, the suffixes c and s have the 
same meaning as before, and L denotes the operator ὃ, 0?/dy” +b,d/dy +b,. 
In deriving the above results it has been assumed that u, Ou/OX, Me ὧν are 
continuous in x and satisfy the conditions of Theorem 8. It is also assumed 
_ that differentiation with respect to y may be interchanged with the integra- 
tion operation in equations (3.7) and (3.8) giving 
a,[—B?U,— (@u/éx), _ ,]+a; U.+LU, =F, (3.9) 
a,[—B?U,+Bu), _,]+a, U,+LU, = F.. (3.10) 
Thus, if du/éx is known at x = 0, equation (3.9) becomes a differential 
equation for U, regarded as a function of y; similarly, if u is prescribed 
atx = 0, equation (3.10) gives the differential equation satisfied by U,. 
Once either U. or U, has been found then Theorem 10 can be used to 
find u. Various assumptions have had to be made about the unknown 
function and it has to be verified that the function eventually determined 
satisfies these conditions. The above general procedure is formally 
identical with that described in ὃ 2.3 for the solution of inhomogeneous 
boundary value problems using Fourier series. 
Ifx ranges from — οὐ to 00 and the conditions onuareu > Oasx > +00, 
then the exponential Fourier transform has to be used. In this case it is 
possible to consider the more general equation 


07 u 


du 
αι 52 14,5 ΟῚ ΜΈΓ τα f, (3.11) 


τω οι) ἐ ον Ty eee οι νόον. 


where α,:. 2. ἃς are constant. Taking the Fourier transform of equation _ 


(3.11) gives, using the result of Problem 3.3, 
—f*a, U+ifBa, U+a,U+LU =F. (3.12) 


Various assumptions concerning continuity and interchanging operators 
have again been made in deriving equation (3.12). This is a differential 
equation for U and, after it has been solved, u may be found from Theorem 
8. 


Problem 3.9 Solve 07u/éx* = du/ét, -- οὐ < x < οὐ, t > Ὁ under the 
conditions u — 0 as |x| -- oo, u(x,0) = exp(— 


Solution. The appropriate transform to use is clearly the exponential 
one, and taking the transform of the equation gives —?(U) = 
δύ δι. At τ = 0, u = exp(—x”) and hence U = [ἢ exp(—iBx— x?) dx 
att = 0,and it follows from Appendix 2, equation (5) that U = πξεχρί-- 18?) 
when t = 0. Solving the differential equation for U gives 


U = A(B)exp(— 71) 
where A can be a function of β. The condition on t = 0 shows that A() 
is equal to the value of U at t = 0 and hence 


U = ntexp[ —487(1+41)]. 


The inverse of x*exp( —+?”) is exp(— — x”) [Appendix 2, (5)] and from Appen- 
dix 2, (1) it follows that the inverse of exp(—+A7?) is A~ ‘exp(—A~?x?). 
Hence setting 2 = (1+4t)* shows that u = (1+4t)~ *exp(— x?)\/(1 +41). 
It can be verified in this case, by direct substitution, that u is the required 
solution. oO 


Problem 3.10 Solve 07u/dx?+07u/dy? =0, O0< x< wm, O< y<d, 
under the conditions u(x,0) = u(x,d)=0, x >0; u(0,y) = 1—y/d, 
0<y<d;u-Oasx > οὐ. 


Solution. This problem could be solved by either representing u as a 
Fourier sine series in y or taking the Fourier sine transform with respect to 
x, and we choose the latter approach. Taking the Fourier sine transform 
of the above equation gives 


os — B(1—y/d), [c.f. equation (3.10)]. 


The general solution of this equation is 
U, = Acosh By+B sinh By +(d—y)/Bd, 
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and the conditions U, = 0 for y= 0 and y =d [which follow from 


u = Oony = 0, 47 yield 
y — sinh BQ-@) (4--ν) 
᾽ β sinh Bd Bd ᾿ 
_ 2 |* sinh B(y—d)sin BxdB 2(d—y) [@ sin Bxdp 
Hence 4) ἜΤ Ganado π᾿ B . 
Setting « = y—d, y = d in Appendix 2, equation (8) gives 


* sinh x(y —d)sin Bx is is Ἄτην" πβ 
| a + " dx = tan tan =o 4 ftann( 3) 


and interchanging the variables B and x gives the first term in u. The 
second term in u can be evaluated by using the identity derived at the end 
of Problem 3.6 (replacing 8 by x) and thus finally 
we few tant con ee 
“= 1 7 ee cot $F tanh | 
Clearly this function satisfies all the conditions imposed. oO 


Problem 3.11 Solve the equation of the previous problem in 0 < y<d, 
x > 0, subject to the conditions u(x,0) = 0, u(x,d) = 1; x > 0; 
u(0,y) =0,0< y<d. 


Solution. The sine transform is still applicable and now 
0°U, 
dy? 
As u vanishes on y = 0, U, must also vanish on y = 0; u(x,d) = 1 and 
hence U, on y = d is the sine transform of unity. However, unity does not 
satisfy the conditions of Theorem 10 and it does not seem possible to 
calculate its Fourier sine transform. It follows however from the last 
identity of Problem 3.6 by replacing β by fx, (x > 0), that there exists an 
F(f) (in this case 2/nB) such that [ὦ [F(f)sin Bx]dB = 1. F(B) is not abso- 
lutely integrable and thus Theorem 10 cannot be used to invert the above 
Tesult. We therefore seek an integral representation for u of the form 
u = [σ᾽ G(B, y)sin Bx dB. Substituting this in the equation for u shows that 
G satisfies the same equation as U, above. The boundary conditions give 
_ G(B,d) = 2/nxB and G(B,o) = 0, and hence G = 2sinh By/B sinh Bd and 
ΤΕ :" sinh By sin Bx dp 
ajo  BsinhBd — 
This integral is of the same form as that of Problem 3.10 with y—d now 
replaced by y and hence 
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—f?U, = 0. 


(For all problems with homogeneous boundary conditions the above 
method of assuming an appropriate integral representation of the solution 
can be used instead of taking the transform of the equation. The two 
methods are entirely equivalent but, provided the relevant transforms 
exist, taking the transform of the equation has the advantage of also being 
applicable for non-homogeneous boundary conditions. This is also the 
method used for solving ordinary differential equations using Laplace 
transforms.) 


Problem 3.12 Solve the equation of Problem 3.9 for x > 0, t>0 
subject to the conditions u ~ Ὁ as x > οὐ, u(x,0) = 0 and éu/éx = 1, 
x = 0,t > 0. 

Solution. In this case the only possibility is the cosine transform and 
—f?U,—1 = 6U,/ét [c.f. equation (3.9)]. The general solution is 

1 
U. = ager Aiere-63 

and, as U. = 0 for t = 0, it is necessary that A = B~*. Thus applying 
Theorem 10 gives 


‘ dp. 


_ 2 [“ cos Bx(exp(—B?t)— 1) 
=e Sree aN 
Appendix 2, equation (12) gives (interchanging x and f) 

JS exp(—aB?)cos BxdB = 4(n/x)*exp(—x?/4a) 
and integrating both sides of this with respect to « from « = Otoa=t 
gives 


ἦν b= EO pede mint [ α΄ texp(—x?/4a) da. 


β' 0 

0 

Making the change of variable v? = x*/4a enables the integral on the 
right-hand side to be rewritten as 


- { τ exp(—v?)/v? dv = 2ttexp(—x?/4t)+2x [ἢ exp(—v?)d», 
on integration by parts. 


Hence = — 2ttexp(—x?/40)/n* + x Erfo(x/2t%), 
where Erfcz, the complementary error function, is defined by 
Erfez = 2x7 ὁ |” exp(—t?)dt. oO 
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ΡΩΝ ΩΝ 
Problem 3.13 Solve 
δι 1éu é?u 
ait a ta 0, 0<r<a,z>0 


under the conditions u(0, z) finite, u(a,z) = Ὁ and u(r, 0) = 


Solution. This problem can be solved by taking the sine transform with 
respect to z and this gives 
eu 1 8 τον 
ér? ars δὲ Fi age 
A particular integral of the equation is 1/8 and the corresponding homo- 
geneous equation is Bessel’s equation of order zero in the variable ifr. 
The appropriate finite solution of this is J,(ifr) which is normally written 


as I(r). 
The solution satisfying the condition on r = a is thus 


* BL 1,(Ba) 
_2;,°) ,  1,(6r) |sin Bz 
Hence “= =| [ ar β ap. 


_ The first term is equal to unity (c.f. last identity of Problem 3.6) and the 
second integral can be transformed by use of Cauchy’s theorem to yield 
the form of solution obtained by using a Fourier—Bessel expansion 
(cf. Chapter 5). From the numerical point of view the infinite integral 
form is as convenient as the infinite series form. oO 


Problem 3.14 Solve 67u/@x? + é?u/éy? = 0, y > 0. all x, under the con- 
ditions u(x,0) = f(x), u+0 as x?+y? - co, where f is a function 
satisfying the conditions of Theorem 8. 


Solution. The two possible methods of solution are the exponential 
Fourier transform with respect to x or a sine transform with respect to y. 
The latter will lead to an inhomogeneous equation for the Fourier trans- 
form and this generally presents more technical difficulties. We therefore 
take the exponential Fourier transform with respect to x and have 
6?U/dy?—B?U = 0. In order to satisfy the conditions at infinity 
U = exp(—|fly) and to satisfy the condition on y = 0, Α(β) must be 
equal to F(A), the Fourier transform of f. Hence U = F(f)exp(—|fly). 

Theorem 8 now gives a formal expression for u but a simpler one is 
obtained by applying Theorem 9 (the convolution theorem) and using 
the result that the inverse of exp(—| β] y) is y/n(x? + y”) [Appendix 2, (4)]. 
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; 4 lie 
~ Ἄν: ψ 
Hence 
ΕΠ σῦς 
-2|" 2X. fe 
Problem 3.15 Solve 
O*u δε 8 
= (), > 0, all x, 
ax* Fire oe me 


under the conditions du/éy = 0, y = 0; (éu/éx),_. = f(x), where f(x) 
is a function satisfying the conditions of Theorem 8. Also u — 0 as x, 
y > 0. 


Solution. This equation is not of the general type that we have considered, 
but by analogy with our previous examples it would seem reasonable to | 
attempt a Fourier transform and use of the exponential form avoids any 
complications with inhomogeneous terms, Thus 

ou 
U —2p* —~ = 0), 
p* Par T oy 
the general form of solution wana to zero as y > © IS 
U =[A+By]exp(—|B|y), the condition on du/éy on y =0 gives 
|B| A = B. The condition on éu/éx on y = 0 gives iBU = F, where F is 
the Fourier transform of f and thus | 


i 
U- gil +1Bl yexp(—|B|y)FB) 
The inverse of the term multiplying F() is 


-:] oP expt |B\y)exp(iBx) dB 


ε: an SE eae inhigh 4. |” sin Bx exp(— By) dB. 
0 B π 0 


π 
From Appendix 2(9) the first term of the above equation is equal to 
(1/x)tan~ ‘(x/y) and Appendix 2(12) shows that the second integral is 
equal to x/n(x?+ 7). The solution for u can now be written down by 
means of theorem 9. Oo 


Problem 3.16 Solve Problem 3.12 when the condition on x = 0 is 
Gu/Ox = hu where h is a constant and u(x.0) = exp(— 


Solution. If h were zero then the appropriate transform would be a 
Fourier cosine transform and for ἢ - coo the Fourier sine transform 
would be appropriate but neither is suitable for other values of h. It 
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r 


me 


Radin. seems worth considering the possibility of using a combination 
of sine and cosine transforms and we consider U* defined by 


‘Bina s 


εὖ." δ απ a dy usin Bx dx+b \- ucos Bx dx, 


' 


_ where a and ῥ᾽ are functions of £. To obtain the differential equation 
satisfied by U* it is co to calculate 


au (asin Bx +bos px) de = ΟΣ ΙΝ 
Ξ 


< 
ty pt 


- Τ. (βα οο8 Bx — Bb sin py dx 


me (asin px+be0 pry (Ba cos Bx — Bb sin pm | > 
=0 


this means the integrated τῳ vanish as x > oo and 


ἃ (asin Bx+bcos py dx = ΕΞ ταβι} — B?U*, 
=0 

and, in view of the boundary condition on u at x = 0, the first term 
vanishes if af = hb and the simplest choice of a and b gives a = ἢ: 
b= B. Multiplying the equation for u by h sin Bx + B cos Bx and integrating 
with respect to x from 0 to oo formally gives 0U*/dt = — f?U*. 

Att = 0, U is exp(—x) and hence 


U* = ᾿ (h sin Bx +B cos βχ)εχρί-- χ) dx = βίμι- ΠΧ - β2), 
_ Thus U* = (1+h)Bexp(—f7t)/(1 +B”). It is now necessary to determine 
_ the inversion formula for this transform, we have 


UF: = ‘a hu sin Bx dx + | usin Bx) dx = | (mu — sin Bx dx, 
0 0 
(on integrating the second term by parts). From this expression it is seen 
that U* is the sine transform of hu—@u/dx and hence 


Ox 


This equation can be rewritten as 


ὃ 2 2 
— 5 (ue) = cet | U*sin Bx dB 


0 


μι--ν - 3" U* sin Bx dB. 
T Jo 


é ͵ . 1: _— ae. 5 ᾿ ᾿ Ε = 4 4 
te ‘a oo . pe! 4 
᾿ - .- : 
» id 4 f , ' . ᾿ Te } +* 
μὰ ᾿ ΟΝ. my yh >! 
3 ‘ ΡῚ is. 


— p? ᾿- (asin Bx +b οο5 Bx)u dx. 
0 


For any of the integrals to converge it is necessary that u > 0 as x > co: 


é “Δι, a ᾿ 7 ΨῃΨ} i « ; t am 
Le is wet . i ey 


᾿ς “ΝΣ ἃ * 
a i a 
ἱ — a, - . . γ᾽ “ F 


and replacing’ x by a variable w and integrating both sides wid feet ‘a Ἂ 
to w from w = x to w = οὐ gives > 


Aj (ἢ sin Bx + B cos Bx)U* F + 

= A τ h2 ey p? B. ; Α 
The above derivation of the inversion formula has made various assump- 
tions concerning interchanging orders of integration, but it can be 
verified that the inversion formula is valid for all functions f(x) satisfying _ Ἵ 
the conditions of Theorem 8. Thus the solution of the present problem ὃ 
can be obtained in the form of an infinite integral. rs ie 


EXERCISES 
1. If f(x) is exp(2x) for x < 0 and zero for x > 0 show that 


fix) = ae 70s is psn dp. 


2. Obtain the Fourier sine transform of exp(—x)sin x, (x > 0) and hence 
show that, for integer n,|* (x sin nzx)/(x* +4) dx = 0. μ᾿. 
3. Solve the integral equation 
=0, 0< Bl, 
j, fedsin Bxdx =1, 1<fB <2, 
=O, 6 .> 2. 


4. Solve δι δι = 07u/éx?, x, t > 0, subject to the conditions 
u(0,t) = sin@t,u > Oas x > c,t > 0; u(x,0) = 0.x > 0. 


5. Solve 6?u/@x? —3u = éu/ét, x,t > Ὁ under the conditions u(x,0) = 
x > 0; du/ox = εἰ χε 0,t > 0; lim u(x,t) = 0,t > 0. 


δ δ: 
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Chapter 4 


Generalized Series Expansions 


4,1 Introduction In Chapter | it was shown that a wide class of functions 
could be represented by series of trigonometric functions and we now 
examine the extension of this principle to represent a function as an 
infinite series of certain other types of functions. In practice it is only 
useful to consider a representation of f on [a,b] by a series of the form 


9) 
>> a,¢, where the φ, satisfy a relation of the form 
tile | 


ρφ͵ φ͵,, ax = 0, πὶ # ἢ, (4.1) 


where p is a given function. In this case the functions ᾧ, are said to be 
orthogonal with weight p on [a,b]. All the Fourier series previously 
considered are of this form with p = 1 (c.f. Problem 1.3), 

If it is assumed that, for a given set of functions , , a function f is such 
that 


f= Σ δ, ἀπ τ ἢ (4.2) 


then, multiplying both sides of this equation by pd, and integrating with 
respect to x from x = a to x = b gives 


Γ᾿ ef b, dx = y γ᾽ ρφ, Φ, ἀν, 


(assuming that the order of summation and integration may be inter- 
changed). Equation (4.1) now shows that 
| [ pf, ax =a, [ pd? dx, (4.3) 


thus yielding an expression for a,. The above analysis is purely formal as 
it need not be true that, for a given set of ¢,, equation (4.2) should hold, 
for any f. (If the @, were chosen to be 1, cos 2x,cos3x,...on [0,2] then 
for f = cos x, equation (4.3) would give a,,=0, thus yielding an obvious 
contradiction and showing that the above set of functions would not be 
suitable to use in representing a function on [0, 2].) In those cases where 
equation (4.2) is meaningful the coefficients defined by equation (4.3) are 
known as the generalised Fourier coefficients of f and the resulting series 
representation as a generalized Fourier representation of f. The formal 
a given in Problem 1.3 is a special case of the method outlined 
above. 
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The above example shows that a necessary condition for a representation 
of the form of the right-hand side of equation (4.3) to be valid is that there 
exists no function orthogonal to all the ¢,. In making statements of this 
nature concerning the existence of a function, one has to be somewhat more 
precise and specify the class of functions considered (e.g. functions con- 
tinuous on an interval, or the class of integrable or square integrable 
functions over an interval). The class of functions is sometimes said to 
constitute a function space. If the functions #, are orthogonal with weight 
p on [a,b] and such that there is no function f of the same class, ortho- 
gonal to all the @, and such that [ἢ pf? dx # 0, then the set of functions 
φ͵ is said to be complete within the given class. If the series of equation 
(4.2) with a, defined by equation (4.3), converges to f, except at possibly 
a finite number of points in [a,b], for all functions f of the class con- 
sidered, then the set @, is said to be closed in the sense of pointwise 
convergence. This terminology is introduced here to make the reader 
aware of it but the concepts will not be used in the text. 

We now examine briefly two general classes of problems which lead 
naturally to the idea of expanding a function in terms of two particular 
sets of functions. We consider first the problem of solving Laplace’s 
equation in a spherical region when axi-symmetric boundary conditions 
are prescribed on one, or two, spherical surfaces. Laplace’s equation 
becomes, on assuming axial symmetry, 


1 ὃ Ou δ." ὦ . ,OUu 
ἀξ» 3 δ 000 (sino) = οι 
where r, θ are the usual spherical polar coordinates. Writing u as Ε (r)F ,(@) 
and using the procedure of separation of variables gives 


rdr\ d 
Wf. dF, 


and, on writing x as cos 6, equation (4.4) becomes 


d ΘΕ, 

ial 5 ee) ee = 0. 4. 

=| x) | + aF, 0 (4.5) 
In physical problems F,, has to be finite for all @ (i.e. |x| < 1) and it can be 
shown from the theory of ordinary differential equations that equation 
(4.5) only has such finite solutions when 2 = n(n+ 1), where n is an integer. 
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In this case there is a finite solution, the Legendre polynomial, P(x), of 
degree n satisfying 


Lt 
| dx 
and the condition P (1) = 1. 
When 4 = n(n+ 1) independent solutions for F, are γ" and r~"~*, hence 


τ os] + n(n+1)P, = 0, (4.6) 
dx 


‘an appropriate series representation for u is νὰ (Α,,"- Β,γ 5 ἪΡ (ςο5 θ). 


Thus if u is prescribed for r = a and r = b in terms of θ and if any function 
of @ can be expressed as a series of Legendre polynomials then A, and 
B, can be found. Boundary value problems of this type lead to the concept 
of expanding a function as a series of Legendre polynomials. The formal 
properties of such an expansion are given in a subsequent section but it is 
convenient to summarise some of the simpler properties of these functions 
at this stage. 
An explicit form for P,(x) is provided by Rodrigues’ formula 


P(x) = 4502-17, (4.7) 

some particular cases of this being 
P(x) = 1, (4.8) 
P(x) =x, (4.9) 
2P,(x) = 3x? —1. (4.10) 

Also 

P(—x) - (- ΤΡ (x), (4.11) 
P, (0) = er n = 1,2,..., (4.12) 
P,_, (0) = ‘ (4.13) 
P,,(0) = (4.14) 
P,,.,(0) = a 1)P,,,(0), (4.15) 
P(0) = Ρ΄(0) =1, (4.16) 
(1—2xt 2). = YP) (4.17) 


The Legendre polynomials are orthogonal with unit weight on [—1, 1] 
and 


: 2 
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_ Laplace’s equation in cylindrical polar coordinates (ρ, Φ, 2) becomes, 
assuming axial symmetry, 


Lo (ou), au _ 
pop Pap dz? 
writing u = G,(p)G,(z) and using separation of variables gives 
ld/ dG, 
+o (oS) 4 AG, = Ὁ, (4.19) 
40 
7 2_ 1G, =0 
Equation (4.19) can be rewritten as 
d*G, 1dG, 
dy? + fi G, = VU, 


where y = A*p; this is Bessel’s equation of order zero with solutions _ 


J,(y) and Y,(y). For finite solutions in boundary value problems for — 
regions including the origin Y, must be excluded. For the particular 
class of problem when u is required to vanish on p = 1 it is necessary 
that J,(4*) = 0 and thus A* must be one of the zeros of J,(y). There are 
an infinite number of these which we shall denote by j,,. The equation for 
G,issoluble in simple terms and thus the appropriate series representation 


is Σ᾽ (A, coshj,, +B, sinhj,, z)Jo(jp,?). If u is prescribed on z = a 
n=0 


and z = b and if any function of p can be expanded as a series of J,(j,,p) 
then A, and B, may be found. 

In § 3.3 the general concept of expanding a function in [0,1] asa series — 
of terms of the form J,(j,,x), where J,(j,,,) = 0, is examined in detail. Some 
of the simpler properties of Bessel functions are summarised below. 

The Bessel function J,(Ax) satisfies the equation 


aj, idJ, y? 
de τα a+(2-B) 55 sia 
_ & (ey 
Also d[ x" J (Ax) ]/dx = Ax’J,_ (Ax), (4.22) 
bo (Ax)]/dx = i (A), (4.23) 
J+J, , (x) = “7,00 (4.24) 
J? se —J.,,. (0) = 25 (x), (4.25) 


where the dash denotes the derivative with respect to the argument. The 
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set of functions J,(j,,x) [n = 1,2,...] are orthogonal with weight x on 
[0, 1] and 


ας x)dx = 437, (ὦ (4.26) 


Series involving J,( j,, x) are generally known as Fourier—Bessel series. 


4.2 Expansion in Terms of Legendre Functions The basic result is given 
by the following theorem. 


Theorem 11 If f is piecewise smooth in [—1,1] then the series 
2. a, P,(x),where 2a, = (2n + 1) Εν f (x)P,(x) dx, converges for |x| « 1 


to a (x) at all points at which f is continuous and to 4{ f(x+0)+f(x— —0)] 
at points of discontinuity. At the end points x = +1, x = —1 the series 
converges to f(1—0) and f(—1+0) respectively. The particular form 
occurring for a, can be deduced directly from equations (4.3) and (4.18). 


1 I n 
| Ff (x)P,(x) dx = τ ᾿" 7 0) (x*—1)"dx, [from equation 4.7] 
-ι 


irr aoe 
~ πη! LL rata e> | - ΕΝ τ στ 
(x* -- 1)" dx, [on integration by parts. ] 
(x? -- 1)" and its first (n—1) derivatives vanish at x = 1,x = —1 and thus 
the integrated part in the above expression vanishes. The process can be 
repeated n a giving 


— 


τ οὐ --1}" dx. (4.27) 


Equation 4.27) is often als in pias coefficients in Legendre 
expansions. 


Problem 4.1 Find the Legendre expansion on [ — 1,1] of x?. 


Solution. Equation (4.27) and Theorem 11 show that the coefficient 
of P (x) is given by 


_ (=1)"2n+1) 


a Batty 


a” (:2V 2 

7.» (ὰ γχχἪ -- 1)" dx. 

The nth derivative of x* is zero for n 2 hence a, = 0, n > 2 and 
1 

α; εὖ [ (x?-1P dx τ, a, = -3 [ χα’ -- ἢ dx = 0, 

᾿ i 

a, = [i x? dx =< 
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Hence x? = 2P,+4P,; this particular result could have been obtained 
more directly from equations (4.8) and (4.10). 


It can be shown similarly that the coefficients of P, in the expansion of 


x” vanish for n > m. In this case the Legendre expansion would only have 
a finite number of terms and is best obtained by comparing the coefficients 


of powers of x in the expansion x” = P (x), where the explicit 


expressions for the Legendre polynomials are substituted on the right- 
hand side. 


Problem 4.2 Expand cos 2¢ as a series in P (cos@?) n0 < @ < π. 


Solution. In this example x has been replaced by cos @ and the same 
substitution should be made in the integral for a, . It is, however, simpler 
in this case to use cos 2¢ = 2cos*¢—1 and equations (4.8) and (4.10). 
Thus 


cos 2p = 2cos*¢—1 = 2[2P,(cos @)+1]—1 
= ΞΡ, (το φ) --ἰ = $P,(cos¢)—4}P,(cos@) (since P, = 1). O 
Problem 4.3 Find the Legendre expansion in [ —1, 1] of f(x): 
f(x) = 9, (-l1 <x <0), f(x) = 1, O<x< 1). 
Solution. In this case a, = 4(2n+1)|} P(x) dx. The integral can be 
evaluated by using equation (4.7) but a more direct approach is to use 
equation (4.6), giving 
(2n+1) », aP, n+ 
dx 


em ~ 2n(n+1) 


= ——— P’(0). 
2n(n+ 1) 0) 


P’ (0) can be calculated from equations (4.13) to (4.16), which show that it 
is zero for n even (n τέ 0) and 


Oy Gd | dk Ce 
ἄγκεα = (4k 422k +2) 2.4.6...2k 
a 


Kite 1,2,.... 


Nil 
, 


i | — 
; as ay = 


By theorem 11 the series should converge to 1 at x = 0, which is a point of 
discontinuity and the above results confirm this as, from equation (4.13), 
P,,,.,(0) = 0 and the value of the series at x = Ois a, = 4. 0 


Problem 4.4 Show that it is possible to represent a function f(x) de- 
fined in [0, 1] as a series of Legendre polynomials P (x) of even order. 


Solution. Since Theorem 11 refers to the interval [ — 1, 1] it is necessary 
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to find a function F(x) equal to f(x) in [0, 1] and such that the coefficients 


of odd order in the Legendre expansion of F are zero, i.e. 


ft F(x)P,(x) dx [; P (x)F(x) ἀχ- \" _ PalX)F(x) dx 


᾿ P (x)LF(x)+(—1)"F(—x)] x dx =0, nodd 


on replacing x by —x in the second integral and using equation (4.11). 


Thus F(—x) = F(x) and hence F(x) has to be defined for — 1 < x < 0 
as f(—x). The coefficients a,, of P(x) are given by 
4n+1 [᾿ - 
a,, = | F(x)P,,(x) dx = 4n+1 ἢ f (x)P,,(x) dx. 
-1 
The series converges at x = 0 to 2LF(0+0)+ F(0— —0)] = f(0). Oo 


Similarly f(x) can be represented as a series of Legendre polynomials 
of odd order by defining F(x) = — f(—x) for -1 < x < 0, 

In this case the resulting series converges at x = 0 to 
AT F(0+0)+ F(O—0)] = 0. 


_ Problem 4.5 Find the series representing x on [0, 1] in terms of Legendre 
polynomials of even order. 


Solution. This is a special case of the previous example with f(x) = x 
and hence a,, = (4n+1) [Ὁ x P,,(x) dx. The integral is best calculated by 
using equation (4.7) giving 


ig 4n+1 ms 
a 2?"(2n)! J , ΩΣ “Juan 


(4n+ 1) it "ως 1" d2"- :χ — ee yt 
= 222m)! )| © axe =f τ Ξε τ κ᾿ 
n # 0. 


The integrated part vanishes at both limits and thus, after a further inte- 
gration, 
_ | 4@a+1)| 4552 ee 
a,, = 2?"(2n)! qyaen a (x a) ᾿ Mg n#0, 
the expression in square brackets vanishes at x = 1 and can be evaluated 


at x = 0 by expanding (x* — 1)" by means of the binomial theorem, yield- 
ing 


(x? — 1)?" dx 


SS ae aS erm 
an" 2*(Qn)!dx™-? & ri(Qn—n! 


x = 0. 
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It is only the term x?"~* that produces a non-zero contribution after 
differentiating 2n—2 times and setting x = 0, and thus 
_ (4n+1(2n—2)(— 1"? 
“an τ τ πῃ — ΝΠ  ἢΙ 
a, can be obtained directly from the initial expression (since P, = 1) 
and is equal to 4. 


n # 0. 


Problem 4.6 Expand (5—4x)~* as a series of Legendre polynomials. 


Solution. The straightforward approach using Theorem 11 will lead 
to evaluating integrals of the form {* (5—4x)~+P_(x) dx and this does not 


seem to be a very simple task. It is, however, worth checking in the case 


of an expression of the form (a+bx)~* whether it is a particular case of 
the left-hand side of equation (4.17). In the present example 


(5-- 4χ) 3 = (4—4x41)"# = 41 —24x 44)-4. 
This suggests applying equation (4.17) with t = 4, giving 
(5—4x)-+ =1 Y 2-"P (x), oO 
n=0 
4.3 Fourier-Bessel Series In the present section we examine the possi- 


bility of expanding a function on [0, 1] as a series of terms proportional to 
J j,,x) where J,(j,,) = 0, and the relevant expansion theorem is: 


Theorem 12 If f(x) is piecewise smooth in [0, 1] then the series 


᾿ a, J (jy_X), where 


n= 1 


a= m5). xf (x)J (ij, x)dx, 


converges to f(x) at all points at which f is continuous and to 

+L f(x+0)+f(x—0)] at points of discontinuity. The series converges 
to zero at x = 1 for all v and for v > 0 converges to zero at x = 0. The 
particular form for a, may be deduced from equations (4.3) and (4.26). 


Problem 4.7 Obtain the representation in [0,1] of unity as a series of 
functions J,(j, x) where J,(j,) = 0. 


Solution. This is a particular case of Theorem 12 with v = 0 and f = 1 
and thus 


pe δς 
a, = Fu) | xJ,(j,x) dx. 
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The integral can be evaluated by applying equation (4.22) with v = 0 
and / = j,, giving 
iz 2 

————— xJ Rs ΕΣ Ἐπ σου [] 
© = TPG) Cee 0 iJ Gy) 
Problem 4.8 If μι. are the roots of J,(4u) = 0, obtain an expansion in a 
series of J,(u,x) of the function f(x) defined on [0,4] by f(x) = 1, 
O<x<2, f(x)=0, 2<x<4 [(2) -- 2. 


Solution. Defining a new variable y by +x gives an expansion problem in 
[0, 1] in the new variable. The j,, of Theorem 12 have now to be replaced 
by 4u, and thus the coefficient a, of Jo(jo, y)L= Jo(u,x)] is 


a 


J,(2u,) 
4u, Ji(4u,) 
(equation (4.22), with v = 1,4 = 4y,). ΓΙ 


2 3 2 a 
Pan) [ yf (4y)J,(4u,,y) dy Flan) |. yJy(4u, νὴ dy 


The only cases where the coefficients can be evaluated are those when 
the integrals occurring are such that the equations (4.22) or (4.23) can be 
used. There are some other special cases where the coefficients take on a 
comparatively simple form, but these cases require a detailed knowledge 
of the properties of Bessel functions. Further expansion problems for 
Fourier—Bessel series and their generalisations occur in Problems 5-2, 
5.4, 5-5, 5-6, 5:9 and 5-10. 


4.4 Sturm-Liouville Theory It is effectively shown in Problem 2.1 that 
the set of functions sin nzx/c(n = 1,2,...) is generated by determining / 
such that X” = —AX has non-trivial solutions in [0, c] vanishing at the 
end points of the interval. This type of problem is known as an eigen 
value problem, the appropriate values of Ἢ are the eigen values and the 
corresponding solutions are known as the eigen functions. Another 
example is provided in Problem 2.3, where it is shown that the functions 
sin (2n—1)rx/2I are eigen functions for the above equation in [0,/] when 
the conditions at the end points are X(0) = X’(I) = 0; the corresponding 
eigen values being (2n — 1)?2?/4/?. 

Both the above examples are special cases of a general class of eigen 
value problem, the Sturm-—Liouville problem, which is such that an 
arbitrary function can be represented in terms of the solutions to the eigen 
value problem. 

The Sturm-Liouville problem in its most general form is the determina- 
tion of the eigen values and the corresponding eigen functions in [a,b] of 
the differential equation: 
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d/{ dd 
τῷ +40 +Apod = 0, (4.28) 


where p, 4. p are given functions and the boundary conditions at the end 
points are either 


d 
1,0+0,° = 0, x = α, (4.29) 


d 
Bot B= 0, x=b (430) 
or 


d d 
k,$(a)+k,9(b) = 0 = m(<) ; +m) aa 


The first class of boundary conditions are known as unmixed boundary 
conditions and the second class are normally called periodic boundary 
conditions. If τέ 0 in [a,b] then the problem is said to be regular. 

The principal results of Sturm—Liouville theory may be summarized in 
the following theorem. 


Theorem 13 Let p, p’, 4ρ and (pp)” be continuous real functions of x on 
[a, δ] with p > 0, p > 0 and the constants «,, «,, B,, B,, Κι, k,, m,, m, 
be real. Then the Sturm—Liouville problem has an infinite number of eigen 
values 4, all real and not more than a finite number being negative. The 
eigen functions ¢,, @,, corresponding to different eigen values ἃ... 4, 
are orthogonal with weight p. 

Any continuous function f(x), with a piecewise continuous first deriva- 
tive, can be represented in a < x < b as an infinite series of the form of 
equation (4.2) with the coefficients defined by equation (4.3). Furthermore 
if the function/ satisfies the same boundary conditions as those of the 
Sturm-Liouville problem then the series converges to f in [a, b]. 

For the case of unmixed boundary conditions (1.6. equations (4.29), 
(4.30)) there cannot be two linear independent eigen functions correspond- 
ing to the same eigen value. Also if the additional conditions x, «, < 0, 
B,B, > 0, are imposed then A, > 0. 

The orthogonality conditions also hold for the unmixed problem 
with p(a) = 0 provided that the equation (4.29) is dropped and with 4 
p(b) = 0 provided that equation (4.30) is dropped. | 

Note. The Legendre and Fourier—Bessel expansions of the previous ) 
sections are particular examples of irregular Sturm—Liouville problems | 
and the orthogonality of these sets of functions is a direct consequence of 
the above statement concerning orthogonality for the particular case when 


| 
ὅδ Ὁ Σ 


πο eye 
"μ᾿ 


p(a) and /or p(b) are zero. The boundary conditions for the Sturm—Liou- 
ville problem for the Legendre expansion are that @ is bounded at x = 1, 
— 1 and those for the Fourier—Bessel expansion are that ᾧ is bounded at 
x = Oand vanishes atx = 1. 

The set sin nzx/c clearly corresponds to the case p = 1 = p,q = Oin 
[0,c] with «, = B, = 0 and the sets cos nzx/c, sin(2n— 1)nx/c, 
cos(2n — 1)mx/2c, are also generated by the same equation with different 
sets of values of a, β [for cosnnx/c, x, = B, = 0; for sin(2n—1)xx/2c, 
a, = B, = 0; for cos(2n—1)xx/2c,a%, = B, = 0]. The various expansions 
quoted in Chapter 1 are particular cases of Theorem 13. 


Problem 4.9 Find the Sturm—Liouville problem in [0,2c] whose eigen 
functions are cos nzx/c, sin ππχίς. 


Solution. These functions satisfy X” +AX = Owith A = n?x?/c?, cos nnx/c, 
sin nmx/c do not satisfy the same mixed conditions at x = 0 and x =c 
but satisfy the periodic conditions ¢(0) = (2c), φ΄0) = ¢'(2c), and we 
thus examine the Sturm-—Liouville problem for the above equation with 
these conditions. For A # 0 the general solution is A cos A*x +B sin J?x 
and the boundary conditions give 


A = Acos2c/*+Bsin 2c/}, B = Bcos2c4*— Asin 2c/}. 


These equations will only have non-zero solutions if their determinant 
vanishes and this gives cos 2cA* = 1 and hence A = n?n?/c?. If this value 
of A is substituted into the equations they degenerate, and hence any values 
of A and B are possible and in particular cos nzx/c, sinnzx/c are both 
eigen functions, illustrating that the result concerning the non-multiplicity 
of eigen functions is not valid for periodic conditions. 

For 4 = 0 the general solution is Cx+D and the conditions give 
D = 2Cc+D,C = C, and hence C = 0 and in this case there is only one 
eigen function, i.e. cos nmx/c (n = 0). [] 


Problem 4.10 Determine the eigen functions and eigen values for 
d’p/dx*+1b = 0,0 < x < 1, with (0) = O and hd(1)+¢'(1) = 0. 


Solution. For 2 τὰ Ὁ the general solution is A cos A*x+BsinZtx and 
as ¢(0) = 0 we have that A = 0; the condition at x = | gives 
hB sin 4* + 2*B cos 4+ = 0. In this case the product f, β, of Theorem 13 
is equal to h and hence, for h > 0, A will be positive thus A? will be real. 
This can also be deduced directly by writing 2 = —u? which gives 
hBsinh u+Bucoshu = 0, and, if B is not identically zero, tanhu = —u/h 
and for h > 0 this is not possible. If h is negative then, as tanhu is a 
monotonic increasing function of u, there may be one solution. The slope 
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of the tanh curve, however, is less than unity and, as both sides of the 
equation are equal for u = 0, it thus follows that the curves will not 
intersect for h > —1. Thus for ἢ < —1 there will be one negative eigen 
value 4 = —u?, where uw is the positive root of tanhu = —u/h and the 
corresponding eigen function is proportional to sinh ux. 

The remaining eigen values are given by 4, = v? where the v, are the 
positive roots of htanv = —v. (It is obvious from sketching the graph of 
tan v that there will be an infinite number of solutions of this equation, 
this is also implied by Theorem 13.) The eigen functions are proportional 
to sin v, x. . 

The general solution for 2 = 0 is ¢ = Cx+D, and since ¢ vanishes at 

= 0 it follows that D must vanish and the condition at x = 1 gives 
Ch+C = 0. Thus for the case h = —1, A = 0 will be an eigen value with 
the corresponding eigen function being proportional to x. O 


Problem 4.11 Determine the eigen values and eigen functions for 
x7d?/ dx? +3xdb/dx+i = 0, in [1,e] with ᾧ vanishing at the end 
points of the interval. 


Solution. This is not in Sturm—Liouville form but can be rewritten in 
that form with p = x°, p = x. The first step is the determination of the 
general solution and the form of the equation suggests considering @ to 
be of the form χ and this gives (u+1)*+A—1 = 0. Thus for ὁ # 1 the 
general solution is 

φ = A(t a8 =1 4 By (bay? = 1, 
The condition at x = 1 gives A = —B and application of the condition 
at x = ὁ gives Asinh(1—/)* = 0. Thus (1—A)* = ἱππ and Δ = 1-- πἦπ2, 
n = 1,2,..., and the eigen functions are proportional to x7 '[x’""—x"""]. 

The above analysis does not show that 2 = 1 is not an eigen value and 
this case has to be treated separately. For A = 1 it can be verified that 
independent solutions of the equation are x~* and x~ "In x, and thus the 
general solution is χ *[C+D1n x]. 7 
The condition at x = 1 gives C = 0, and from ¢(e) = 1 it follows that 
D = 0. Hence A = 1 is not an eigen value. oO 


Problem 4.12 Obtain, for the case h > 0, an expansion of unity in [0, 1] 
in terms of the eigen functions of Problem 4.10. 


Solution. The eigen functions are proportional to sinv,x where 
htanv, = —v, and thus, by Theorem 13 and equation (4.3), 
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where {, sinv, x dx = a, J} sin’v, x dx = 4a, [ (1—cos2v, x) dx, hence 
a, = 4(1—cosv,)/(2—sin 2v ). : oO 


In using equation (4.3) we have assumed the result that the eigen 
functions of a Sturm-—Liouville problem are orthogonal, and in this 
simple case we shall give a direct verification of this general result. 
Integrating by parts, we have | 


; Ln υ 
[sin gesin, x ax = 5 ne cost ἘτῸ COS UV, X COS DV, x dx 
0 


δ᾽ ἢ υ 
= --- — πιο 9 
sin v,, COS v, 2 sin D, COS v,, 
n n 
2 
υ; ; : 
+-} | sinv,xsinv, x dx. 
v, 0 τ 
Hence 
v\ (2 
ck : } COS V_ COS L v 
( ] | sin v, x sin v,, x dx = ——®——*| tanv, ——"tan pv, |. 
n 0 D,, - 5 


As v,, v,, both satisfy (tan v)/v = —1/h the right-hand side vanishes and 
hence, for v,, # v,, so does the integral on the left-hand side. 


EXERCISES 
Expand x° in [ —1, 1] as a series of Legendre polynomials. 
2. Expand as a series of Legendre polynomials f(x) defined by: 
f(x) =0,-l<x<af(x)=1la<x <1. 
[Hint: Use (2n+1)P (x) = P,, (x)—P, _ ,(x).] 
3. Obtain the Fourier—Bessel expansion of x? on [0, 1]. 


4. Obtain the equation for the eigen values of d?y/dx? 
Ἢ y/dx* + dy = 0, subject 
to the conditions (0)+ γ(0) = 0 = 2y(x)+3y'(z). 4 
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Chapter 5 


Solution of Boundary-Value Problems by Means 
of Generalized Fourier Expansions 


It is possible to use the results of Chapter 4 to extend the methods of 
Chapter 2 to problems when the functions a,, a,, a, are not constant and 
when a linear combination of u and 0u/éx is prescribed on x = 0 and 
x = c. Problems where the equation has the form of equation (2.2) and 
a linear combination of u and @u/éx is required to vanish at x = 0 and 
x = c will be termed homogeneous problems and, as in Chapter 2, the 
cases of homogeneous and non-homogeneous problems will be considered 
separately. 


5.1 Homogeneous Problems The technique is a simple generalization 
of that of §2.2 A product solution XY is assumed and separation of 
variables used to determine one of the functions (which will be assumed 
to be X) as the solution of an eigen value problem and the principle of 
superposition is then used to obtain a general representation. If the eigen 
value problem is of a type such that any function can be expanded in 
terms of the eigen functions then the boundary conditions and equation 
for the corresponding Y can be obtained as for Fourier coefficients. It is 
possible to avoid some of the manipulative work of separation of variables 
in special cases by remembering the general form of product solutions of 
certain equations. It is shown in §4.1 that the general form of product 
solution of Laplace’s equation finite for all @ in a spherical region is 
Y (A,r"+B r-"~*)P, (cos 0). 
n=0 

Another useful result which follows from the analysis of §4.1 is that 
[AJ,(4p) + BYo(Ap)][C cosh 2z + D sinh 4z], where Y, is the second solution 
of Bessel’s equation, satisfies the axisymmetric Laplace equation. 


Problem 5.1 Find a solution u of Laplace’s equation, valid outside the 
sphere of radius a, tending to zero at infinity and such that on the sphere 
u = cos@, where @ is the angle between the radius vector and a fixed 
direction. Obtain also the solution to the problem when the condition on 
the sphere is u = (a?+ f*—2af οο5 θ) ἢ, f >a. 

[The first problem corresponds to determining the total electrostatic 
potential when an earthed sphere is placed in a field of uniform intensity. 
The second problem arises when the uniform field is replaced by that of 
a point charge at a distance f from the centre of the sphere. ] 
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Solution. If spherical polar coordinates (r,0) are chosen with origin at _ 


the centre of the sphere, then the analysis of § 4.1 shows that an appropriate 
expansion to choose for u is 


“= Σ, (A.r+Br* ")P. (cos θ). 


This form automatically satisfies Laplace’s equation and is finite for all 
real θ. In order to satisfy the condition at infinity A, = 0 and the condition 
onr = a gives 


ἱν 9] 


Σ B,a~"~'P (cos@) = οοβθ. 
0 


Thus we need to expand cos@ as a series of Legendre polynomials, but 
cos @ = P,(cos θ), hence B, = 0,n τὰ 1, B, = a’. 
To solve the second problem we require the expansion of 
(a* + f?—2af cos @)~* = [1 +(a/ f)? —2(a/f)cos 07 
and equation (4.17) shows that 


[1+(a/f)? -- 2(α [γος 61" = J (a/fP,(cos 6, 


hence B. = a*"*'/f"*!. The solution is thus 


fie a =) a 

u=— Ὁ =| — ) P,(cos0) = —[r? +a*/f? —2(a?r/f)cos 96] ἢ, 

Σ “(5 s(cos 8) = =[r? +a*/ f? —2(a?r/ f)cos 0) 
n=0 

from equation (4.17). The quantity in square brackets is the inverse 


distance from the point (r,@) to (a?/f,0) and the above form of solution 
is that obtained directly in the electrostatic problem by the method of 


images. " 
Problem 5.2 Solve ᾿ 

δω ldu  07y 

Or ra of? O0O<r<at>dO, 


subject to the conditions 

u(0,t) bounded, u(a,t) = 0, ἐ > 0, u(r,0) = (1—r?/a?), Gu(r,t)/at = 0. 
t=0,0<r<a. 

(This corresponds to determining the displacement at any time of a 
circular membrane of radius a which is fixed round its edge, the initial 
displacement is 1 ~r?/a? and the membrane is initially at rest.) 


Solution. We assume a solution of the form F (r)G(t) and the method of 
separation of variables gives 
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The equation for F is Bessel’s equation of order zero, in order that u(0,t) i 4 


be bounded F must be /,(Ar) and the condition on r = a gives (λα) = 0. 
Thus Aa = j,, and to avoid unduly complicated expressions we drop the 
suffix zero. Thus the appropriate form of solution is 


—_ Ἴ, Jj, "Ια ΓΑ,, cos j, t/a+ B, sin j, t/a]. 


Since 6u/ét = 0 for t = 0 it follows that B, = 0 and thus A, is the coeffi- 
cient of J,(j,r/a) in the expansion of 1—r?/a’ in terms of J,(j,r/a) or, 
equivalently, of /,(j,x) in the Fourier—Bessel series in [0,1] for (1—x*) 
(this follows on writing r/a as x). Theorem 12 gives 


1 
A, = = | (1—x?)x J,(j,x) dx. 
iwn 0 
The first part of this integral has already been evaluated in Problem 4.7 
and the second part can be evaluated as follows: 


1 1 : 
| x? J Linx) dx: = £ | δ Ξ [xJ(j,x)] dx, from (4.22), with v = 1, 
Jn Jo 


0 n 
: 1 
= ate τ | x? J,(j,x) ax, 
n n 0 
on integration by parts, 
ie κ Γ 
= ae "- | Pelt J abi x)] dx, 
n n 0 


from (4.22), with v = 2, 
J U,) 2J fj n ) 


a re 


Thus from the above result and that of Problem 4.7, 
τὸ 4..(}) 

᾿ς aU) 
this can be further simplified by using equation (4.24) and the fact that 
J(j,) = 0 and the final result is A, = 8/j? J,(/,). Π 
Problem 5.3 Solve δι δὶ = 6*u/@x*, 0< x <1, 10, under the 
conditions u(0,t) = 0, du/éx +hu = 0,x = 1,t > 0; u(x,0) = 0, 
0 < x < 1,h being a positive constant. 


Solution. The occurrence of 67u/dx? suggests a Fourier series in x but 
71 


wT = “er —— ——— Cl TY! ee ΨΡ.,  ΡΥΝ var Νὰ — “~- 
#* . 


the condition on x = 1 is not appropriate for Fourier series, We therefore 
assume product solutions of the form X(x)T(t) and carrying out the normal 
separation of variables procedure gives 

X" +AX = 0, T'+AT = 0. 
The conditions on x=0 and x = 1 give X(0) = 0, X(1)+h X(1) = 0 
and thus the eigen value problem for X is precisely that of Problem 4.10. 
Hence, using the notation of that problem, 

u= δ᾽ A,exp(—v?t)sinv, x, 

n=1 

and on setting t = 0 it follows that the A, are the generalized Fourier 
coefficients of x. The weight p in this case is unity and thus 


rs) ΕἾ 
A | sin?v.xdx = | xsinv χ ἄχ, 


i.e. v,(2v, —sin 2v,)A, = 4(sin v,—v, cos v,). oO 
Problem 5.4 Solve 

0*u ldu 06*u 

Sat 2, 5. 5. 0 εν θεν <4, 


subject to the conditions u(0, z) finite, u(r,0) = 1, éu(r,z)/ér+yu = 0, 
r = a; Ou(r,z)/dz+ Bu = 0, 2 = h, where β, 7 are positive constants. 

[One physical interpretation of this problem is the determination of 
the steady state temperature in a cylinder when one of its plane faces is 
maintained at unit temperature while a radiation condition is imposed 
on the other surfaces. | 


Solution. The application of separation of variables to this equation 
has been considered in §4.1 and appropriate product solutions satisfying 
the finiteness condition are J,(Ar)[A cosh A(z—h)+ B sinh A(z—h)]. Since 
homogeneous conditions are imposed on z = h the algebra is simplified 
slightly by using hyperbolic functions of z—h. The boundary condition 
onr = a gives 


J,(Ar)+y J,(ar) = 0, r= a4, 
1.6. λ J,(Aa)+y J,(Aa) = 0, 


where the prime denotes derivative with respect to the argument. 
[J,(z)= — J,(z)]. Thus Aa will have to satisfy the equation μ" ol) +6J,(u) = 0 
where ὃ = ya. It can be shown that there are an infinite number of HN 
all positive, of this equation and of its generalization p J/(u)+6 J,(u) = 
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It can also be shown that any function f(x) continuous in [0, 1] can be 
represented by a series of the form Po a, J,(u,,x) where μὲ, are the roots 


of the above equation. The coefficients. ἃ are given by(G. Ν. Watson Theory 
of Bessel Functions, C.U.P., 1952, Chapter 18) 


242, J) FO, iy, ὃ αἱ 
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This more general type of series is known as a Dini series and is an obvious 
generalization of the Fourier—Bessel series (which corresponds to the case 
ὃ — οο.). For the particular case y +d = 0 the correct expansion for a 
function as a Dini series is (G. N. Watson, Theory of Bessel Functions, 
C.U.P., 1952, Chapter 18) 


Ὁ ν i v+ 
YJ +204 De” |) δ  ὙῸ at 


Reverting to the original boundary value problem it is now seen that the 
appropriate representation is 


J(u, r/a)[A, cosh p,(2—h)/a+ B, sinh p,(2—h)/al, 
n=1 


a 


where the suffix zero is dropped to avoid unnecessarily complicated forms. 
In order to satisfy the condition at z = ἢ, u,B,+faA, =0 and the 
condition on z = 0 then gives 


Σ J (u,1/a)A,[ cosh μι, h/a+(Ba/p,)sinh μ᾽, h/a| = 1, 


and A, can now be calculated from the quoted result for Dini series, the 
integral being evaluated as in Problem 4.7. The final result is 


A [u, cosh μ, h/a+ Basinh μ, h/a] = 2J,(u,)/Jo(u,) + J(u): oO 
Problem 5.5 Solve 
δι Lou _ ou 
Or? rér δι᾿ 
under the conditions u(0, ἢ) finite, du(r, ΠΟΥ = 0, r = 1,t > 0; μίν, 0) = r, 
O< r<ih 


t->0.0;< 7 < I, 


Solution. Following the standard separation of variables procedure 
shows that there are separable solutions of the form 


{AJ,(Ar)+ BY,(Ar)}exp —A7t. 
The finiteness condition gives B = 0 and the condition on r = 1 gives 
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ἣ . | ne i ; 
¢ Te ie a by z 
» ὅν. [ 


mo ‘ieee 


Sela) = Ὁ Thus the set of values for Aare ἢ, ΦΌΝΟΥΣ ἘΞ ie newest. 


ρον οἵ J5(x) = 0 and A can also be zero. The general form for u is 
3 Ἵ mee oa a, J o(B, ryexp(— B71) 
and the condition on t = 0 gives 
ἊΝ a+ Ἔ a, J(B,r) =r. 


The expansion problem is thus a particular case of the Dini series of the 
previous problem and it is also one of the special cases when v + 6 = 0. 
Thus 


1 
4 4 2 1 
ΞΞ 2 -- 2 - 2 - 
-. a, a 3 | 3 ‘i ἴω J5(B,) Ι ᾿ TolP ᾿ - _ 


0 


5 In this case the integral for a, cannot be simplified. 
a Problem 5.6 Solve 


a Of Ou 07u 
o(xH) =F. EE gi sh τῷ ἡ ἐ»Ὁ 


under the conditions u(0,t) bounded, μ(],1) = 0, t > 0; u(x,0) = f(x), 

_ where f isa given continuous function of x, du(x, t)/ét = 0,t = 0,0< x <1, 

(This problem arises in the investigation of the small oscillations in the 
vertical plane of a heavy chain displaced from rest). 


Solution. This equation is not one for which the product forms of 
solution can be instantly recognized, but writing u as X(x)T(t) gives 


q df ax d*T 
ss axae)+ AX = 0, Gi tat = 0. 


: X is required to be finite at x = 0 and zero at x = ἰ and the problem of 
determining X is a singular Sturm—Liouville problem. A changé to the 
y A independent variable y(= x*) gives 


ae {7 dx 
14/2) aux = 


@ 


This is Bessel’s equation of order zero in the variable 224y and, reverting 

back to x, the solution finite at x = 0 is J o(24*x*). The condition at x = | 

_ gives 2(Al)* = j, where the jy, are defined as in Problem 5.2, and thus the 
appropriate form of solution is 


“= ) 4, cos + B, sin J oLi,(x/D*]. 


δὰ n=] 


-. ———————— Υ; << a Oy -»- ote os 
' ᾿ ~ 


The condition on éu/dt att = 0 shows that B, = O and the other condition 
att = O gives 


fo= ¥ A, JoLix/0*, 


This equation is not in a form such that: the Fourier—Bessel theorem can 
be applied immediately but this can be remedied by writing (χ ἢ as z 
and it is then seen that the A, are the coefficients in the Fourier—Bessel 
expansion of f(z7/),0 « 2 < lie. 


A, = Pi πη]. zf (Iz? )Jo(inz)dz = mI. ΠΣ a ne) ax oO 


Problem 5.7 Solve (I? —x*)07u/6x? —2xéu/dx = 67u/dt?, 0<x <I, 
t > 0, under the conditions u(0,t) = 0, μί!, ἢ bounded, t > 0; u(x,0) = x3, 
Ou(x, ἢ δὲ = 0, t = 0, 0 < x < Ll. (This problem occurs when investigat- 
ing the small transverse vibrations of a string rotating about a vertical 
axis.) 


Solution. There will exist product solutions of the form X(x) T(t) where 
(I? —x?)X” —2xX'+AX = 0, T"+AT = 0. 

The equation for X is Legendre’s equation in the variable x/I but only 
positive values of the argument are considered. If X* is defined by X(x) 
for x > 0 and by — X(—x) for x < 0 then Χ satisfies the same equation 
as X, is continuous for all x with |x| < |. X* is also bounded at x = +1 
and thus must be a Legendre polynomial of the variable x/l. Further, 
since X* is odd in x it must be a Legendre polynomial of odd order, i.e. 
for x > 0 X* =X =P, _, (x/l) and A = 2n(2n—1), n = 1,.... The 
general form of solution is thus 


= δὴ {A,,-, Cos [2n(2n—1)]#t -Ὁ Β,...., sin[2n(2n—1)]*t}P,,_ ,(x/D, 


the condition on 6u/ét gives B, _ 
it is found that 


, = 0, and from the value of u(x, 0) 


χ᾽ = ΡΣ A,,~1 Ἐχ,.- (χα. 


The A,,_, could be found by an extension of the anise of Problem 4.4 
τ it is easier in this case to use the result 2P 4(x) = 5x? —3x, ie. 
x? = 1[2P,(x)+3P, ee and 
= §l°[2P,(x/l)+ 3P,(x/D)], 
5.4. = 3P/5,A, = oA, = nx 1,2: oO 
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ἐν λίθος οἷ hn aad ἕω ; “a hs Ὁ a κὰ 
—— τ᾽ ᾿ , 
Problem 5.5 Obtain a solution of [ἀρ] δοοβ equation within the sphere 
r < a, finite everywhere within this region, equal to unity on the sphere 
- for 0 < @ < a and zero on the remainder of the sphere (r, 0 being the 
usual spherical polar coordinates). 

Solution. The appropriate form of solution to use is that of Problem 4.1, 
in this case however the finiteness condition within the sphere requires 
that B, = 0 and the appropriate form of solution is δ᾽ A,r"P,(cos @) 
with n=0 

τα | lL, ἢ <= 6 <4, 

Σ Α, α"Ρ (cos @) = ree pee 
The A, can now be obtained from Theorem 11 on writing x as cos 6; 
this gives A, a” = 3(2n+ 1) be P (x) dx. For n = 0 the integral is equal 
to 1—cosa and for ἡ ¥ 0 it can be evaluated by use of the recurrence 
relation 

(2n+1)P,(x) =P, (x)—P’,_ τοῦ), 
where the dash denotes the derivative with respect to the argument. Hence 


A,a" =3[P_,,(cosa)—P,_,(cosa)], n2>1. oO 
Problem 5.9 Solve 
δι 1du du 
----.ο.. ——_—_—_— τῷ --Ξ-. h, b, 
ai tS a ap θ. << Ria < F< 


under the conditions u(a, t) = u(b,t) = 0,t > O;andu(r,0) = la<r<b. 


Solution. Separation of variables shows that there will be product 
solutions of the form {A Jo(Ar)+ BY o(Ar)}exp(—A7t) and the conditions on 
r=aandr = b give 

A Jo(4a)+ BY,(Aa) = 0, A J,(Ab)+ BY,(Ab) = 0. 
For non-zero A and B the determinant must vanish and thus 
J (Aa) Y,(Ab) — J,(Ab) Y,(Aa) = 0, 


which is the equation satisfied by A. 

An equivalent method of calculating 2 would be to use the fact that 
there will be product solutions of the form R(r)exp(—A?t) where 
R"+R’/r+/?R = 0 and R(a) = R(b) = 0. The problem looked at from 
this point of view is of standard Sturm-—Liouville type and from this 
general theory it is deduced that there will be an infinite number of values 
for ἃ (which we shall denote by 6,, n = 1,...,00) and the eigen functions 
will be orthogonal to each other with weight r. The nth eigen function is 
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J,(6,") Y,(6,a) — ¥4(6,r)J,(6,a) = ¢,]. 
The general form for u is thus 


Σ᾽ a, exp(—50)[J (δ, )Y,(5, a)— Υ (δ, r)J (5, a)], 
1 


where the a, are determined by δ᾽ a,¢, = 1. On using the orthogonality 
n= 1 
property of @,, this gives 


a, |’ γφῇ dr = [ r b,(r) ἀν. 


The evaluation of the integrals in the above expression is somewhat 
complicated, but it can be completed by using various relations involving 
Bessel functions [equations (1) and (10) of G. N. Watson, Theory of Bessel 
Functions, C.U.P., 1952, Chapter 5]. The final results are somewhat 
complicated and hence will not be quoted. oO 


5.2 Inhomogeneous Problems The extension οἵ the “nethods of the 3 
previous section to inhomogeneous problems will now be considered. The 


first step is to determine the set of functions ¢, which would be most 
appropriate to use, these can be found from considering the corresponding 
homogeneous problem as illustrated in the previous problems. The 


form of the solution is thus Σ Y (y)@,(x) and the next step is determining 
the equation for Y,. The Snalont will be restricted to equations of the form 
ie 
p(x) Ox 
on using the method of separation of variables it can be verified that the 


problem of solving the homogeneous problem reduces to a Sturm— 
Liouville problem and that p(x) is the weight function. The above equation 


δ 85 
(x) =) +a,(xu+ bi) ΟΣ +b)" + bm = f 


is then multiplied by pd, and integrated with respect to x from x = 0 to 


x = c. Integration by parts twice on the left hand side will give that this 
side of the equation will involve derivatives with respect to y of 

J, eu b,(x) [i.e. ¥,(y) [ς eb, dx ]dx and the prescribed boundary values. The 
right-hand side will be the known function [ς p¢, f(x, y) dx; thusa differen- 
tial equation will be obtained for the Y,. For the case when the boundary 
conditions are homogeneous the equation for the Y can be obtained more 
simply by expanding f in terms of the appropriate ¢. Substitution of this 
expansion and the corresponding expansion for u into the differential 
equation and equating coefficients of @, will then give the differential 
equation for Y.. 
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fn elliptic problems, when conditions are prescribed for two values of 
x and two values of y and one of the two sets of conditions is homo- 
geneous then it is best to use an expansion which satisfies these conditions 
automatically. 

It should again be remembered that the simplest method, if possible, is 
to reduce the inhomogeneous problem to a homogeneous one. 


Problem 5.10 Solve 
67u ldu é7u 
oF rar oF 

under the conditions u = δε δι = 0 att = 0, u = 0, r = a; u finite for 

all r,0 < r < a,aqm not being a zero of J,(z). 


= —q sin wt, 0 ar<iad& t> 0, 


Solution. The boundary conditions are homogeneous and it follows 
(on considering the case q = 0) by direct comparison with Problem 5.2 
that the appropriate expansion is u = δ᾽ Jo(j,r/a)W,(t) where we use 

n=1 


the notation of Problem 5.2. The appropriate expansion for the right- 
hand side of the equation is obtained from Problem 4.7. Substitution in 
the equation and equating coefficients of Jo(j,r/a) gives 


2q sin wot 
Jn JU,) 
The conditions on u, 6u/ét at t = 0 show that y, = wi, = 0, 1 = 0, and, 


provided j? # w7a’, the solution of the above equation satisfying these 
conditions is 


Ἢ 
Un tan = 


2qa* [j, Shee, 


CO 
PIG)  i,-0'a 


Ψψ, = 


Problem 5.11 Solve 67u/dx? = du/dt, O< x <a, t>0, under the 
conditions u(x,0) = 1,0 < x < a; u(a,t) = 1, du/dx = hu, x = 0,t > 0; 
h> 0. 


Solution. The simplest approach is to look for a solution which satisfies 
the conditions at x = 0 and x = a. Such a solution is (1+hx)/(1+ha) 
and writing u as (1+hx)/(1+ha)+v shows that v satisfies the same equa- 
tion as u and v satisfies homogeneous conditions on x = 0 and x = a; 
this problem can then be solved by the previous methods. 

We shall, however, consider the application of the direct technique 
described above. The homogeneous problem is that when the boundary 
condition on x = a is u = 0, and by comparison with Problem 5.3 the 


78 


ἃ 
᾿ ; 


appropriate expansion for the homogeneous problem is | aa 
Ἔ X (x)exp(—A,t) where Χ΄ Ἐλ,Χ = 0, ἀπά X’ = hX,x =0; X = = 0, Ὑ 

x == a. bs 
The solution for X,, satisfying the condition on x = a is B sin A3(x— τῷ Ἢ Ἢ 


and the condition on x = a gives A*cos λῖα = —hsin Ata. Thus Jta = 
where w, are the roots of hatanw = — w. The form of solution to τ: ἐ 
ea 


sought is thusu = 5° ,()sinw,(x—a)/a and 
n=1 


Ψ, ᾿ sin?w (χ --α)α ἀχ = |: u(x, t)sin w, (x —a)/a dx. 


It can be shown by integration by parts twice and using the boundnt a 
condition on x = 0, a and the definition of w, that 


2 0*u . w(x—a) W, we [Ὁ  . w,(x—a) 
[Sain Hg, ra σε [usin ay, 


Hence multiplying both sides of the partial differential equation by — 
sin w,(x—a)/a and integrating from x = 0 to x = a gives, on using the 
definition of w,, in terms of u, 


2 a 
we = -αῖψον, 4} sin?w,(x—a)/a dx 


0 “" 
w? 4w? | 
ears Ψ, ---- 4 . 
a a*(2w, — sin 2w,) 
wt 4 
Th =A eS ἐξ ΝΑ  ΞΞΕΕΣ 
τ ὃ , exr( ss) (2w, —sin 2w,) 


The condition ont = Ogives1 = δ᾽ w(O)sinw,(x—a)/a, ice. 
n=1 


sin ee?) dx = w,(0) inte Sa, 
a : a 


0 


or * (cos ν -ἢ Ξ aw (θ)γῶνν, -- sin 2w,) 

Ww, 4w,, 

4cos w 
Hence A. = ———_——_.. 
"᾿ς (2w,—sin 2w,) a 

Problem 5.12 Solve 

62 Ou Lou, 0*u ἘΠῚ ὃ ὴ ; 

Or ror. oz > a ee 


‘ . a ᾿ 
b Ἢ βι 8 ' 4 
: . “> rr a ὴ ἂν ἡ ! ᾿ ᾿ *) ae - a 
' \ ee , - i ἢ » } : ἂν": 4 ἡ ν᾿ 
ΜΕΝ a δ « “ hy ᾿ ᾿ ω 5 5 Π 
2S - ὶ > J < = m 8 Ψ ᾿ , ᾿ ΒΥ 
,ε, ᾿ ᾿ > 2 ᾿ .Ἰ a e 4 A > 4 ;" 


Seay se F eis = oe re — ΨΩ ~~ rt ae ‘ a nas Soe 
εν phd se tng AH ee Pi eet a ei Sa Tien γῆς ee Ὁ ὧν; τυ ᾿ 1s e's he) ae | 
‘ “under the conditions ἫΝ loli r=a aad πὸ to cos?0 on the ese of the sphere, r, nO being 
the usual polar coordinates. ’ 
(0, 2) finite, = oe O<z<h; ulr, 0) =U, ulr,h) = u, where | Obtain also the solution of Laplace’s equation valid outside the sphere, Ὁ ne 
Ε ὦ, δου 5 AR ΟΘΉΝΙΒΗ ᾿ tending to zero at infinity and taking the same values on the boundary. 


Solution. The boundary conditions on z = 0 and z = hare simple so 


ἃ ‘that it is worth considering the possibility of obtaining homogeneous a ae i : ᾿αἹ ᾿ 
conditions at z = 0 and z = ἢ. The function u,(1—z/h)+u, z/h satisfies | Oru lou a Oru δζγκα εἰν ὃ | fv 
 theconditionsatz = Oandz = h,andifwis written asu,(1—z/h)+ men v or ror δι᾿ ; 5 a c 

_ then v will satisfy the same equation as u; v vanishes at z = 0 and z . 


subject vig the conditions éu/dér = 0, r = a, u(0,t) finite, t > 0; u(r,0) = 


᾿ ieee is equal to u, —u, z/h—u,(1—z/h) onr = a. δι ὃι = r?,t = 0,0 <r <a. 


This problem for v is of the standard type encountered in the previous 
__ chapter and can be solved by assuming a Fourier sine series representation. 
hs An alternative method of solution would be to attempt to obtain a 
series representation for u which would automatically satisfy the condition 
ξ onr = a. One such solution is u, ; writing u as w+u, shows that w satisfies 
_ the same equation as u, vanishes on r = a and is equal to u,—u, and 
Ομ, το on z =0 and z =h respectively. The problem for w can be 
solved by use of separation of variables and by analogy with previous 
_ examples it can be found that the appropriate expansion is the Fourier— 
Ps Bessel series. O 
In examples of this type where two approaches are possible, it is pre- 
Ἵ _ ferable to use the one with least complications. In particular, if one set of 
Ε - conditions are homogeneous one should use the approach which auto- 
, matically satisfies these conditions. Thus if u, and u, had been zero the 
_ sine series approach would be the less complicated, and for u, = 0 use 
of the Fourier—Bessel method would have advantages. There is the 
additional point that integrals involving Bessel functions are usually not 
as simple as those involving trigonometric functions, and this could mean 
_ that even for non-homogeneous conditions it might be better to use 
trigonometric series where possible. 


4. Solve du/at = du/éx, 0 <x < |, t > 0, subject to the conditions — ῳ 
6u/dx = hu, x = 0, du/dx = —hu, x = |, t > 0; u(x,0) = 1, 0 ΧΕ 
h being a constant. ‘i 


5. Solve the boundary value problem of the previous exercise when the 
conditions on x = 0, x -- ἰ are δι οχ = hu—hu,, 0u/éx = —hu+hu, t 
respectively where u, is a constant. 


᾿ " 
iN EXERCISES 


67u ldu du 
ape toa πα =O, 0<r<a,z>Q0, 


subject to the conditions u(0,z) finite, u(a,z) = Ὁ, z > 0, u(r,0) = 1, 
O<r<au-O0asz- oo. 


_ 2. Obtain a solution of Laplace’s equation finite everywhere within the 
χὰ. “- 


Chapte 4 
2. 26/p*+4. 


ἘΠ x — cos 2x)/nx. 


4. : ἘΞ = | sin ot — βὼ cos +f exe 6] dp. 


(The first two terms can be evaluated explicitly using Tables of Fousere . 


| Ἔ 2 4 © (-1)'* !cosnx transforms (Erdelyi, A., et al. Tables of Integral Transforms, McGraw- ¥ 
Ε΄ ΘΟῊΝ = == “9 i Etat 1" O<x<e. | Hill, 1954, vol. 1, p. 66, nos. (23) and (24)) and are equal to ee 
- agp ‘ exp[ — x(o/2)*]sin[ ot — x(@/2)*].) »"- 
4. τ [(—1)"*1(2n—1)rexp(c)—2c] cos(2n— Ix x a 
Ε pS SN ee eS. ee eee eee 2 ) 
- τ (2η -- 1)π2 +4c? 2c 5. » Ὅτ ΤΣ Γοκρί-- Bt) —exp(—0)] ap. a, 
a 0 he 
[2(—1)"* *cexp c+(2n—1)r] sin(2n— ἢπχ. ak 
ς ἢ Ἔ [(2η. -- 1)2π2-- 4.3] 2ς = ΣΙ 
" ᾿ 8 a 3 f 
ς 2. 4 y cos 2nx 1. ΕἸ οὐ ΞΡ: | 
τς ο΄ -ὶ 2. See Problem 5.8. 
Joli, χ) 
ater 2 ae y 27 0. 1} i(j,)— 2/ ,(j,)],using notation of Problem 5.2. 
| a=1 In 
τῷ ρος am a A to saa | eee ies “λυ da cet 4. (2+34)sin λὲπ- J4c0s λὲπ = 0. | 
| 
a=1 με 
Chapter 5 εἰ 
: παν ἢ... ) 2 2n0 Go 
2. psind, 5p sin θ ee y ie (see Problem 1.6). | 1. (The notation of Problem 5.2 is used in 1 and 3.) . 
: ΞῚ bac P " 
4 y sinh[(2n — 1)xy/2a]cosech[(2n — 1)xb/2a]sin[(2n —1)nx/2a] , 2 ᾿ ae exp(—j,2/a). 
nm os (2n—1) Poe 1 "4 
5- 2,-.4 32 
_ 4, u-—1isa solution of exercise 1. | Ae rates Pier enes Se. F feos d+ ar yond, 
| ae “Ὁ ἋΣ 3. 2a δ J i, γα}, ?J ,(d,)cosj, t/a+a°[j, J,(d,) a 
a+(b— a) ~ po sin[(2n — 1)xx/l]exp[ —(2n—1)?2?t/I?](2n—1)7* 0 " πο 
π εὐ —2J,(j,) sin d, t/a]/d? 75(4.). +3 
. 21 3 aie hatha Ων — fA, sind, 1+h—hcos 2,1 ][cosA,x+h4~'sind, x] exp(—2? ἢ 
Ν hate (—1)"* 'sin(nnx/Dexp[| —n*n*t/I*]n~*. 4. rs Tho, a 2 
- a? 81 = (-1)"* ,cos[(2n — 1)xt/21 ]sin[(2n— L)ax/21] where A, are the roots of (A? —h?)tan Al = 2hd. | x 
XP Σ (2η. -- 1) 5. Writingu as ἡ. +v gives that v 15 (1 —u,) times the solution of exercise 4, _ ᾿ 
5 y [nzsin t—Isin nxt/T] sin nnx - 
(n?n? -- |?) Ϊ - 


n=1 


i 
- αἱ 


Ftp) =|" exp(— ips) ἃ = | 


λσαλ F(AB) at ;- τ' 

F(x) 2nf(—B) = 

{(λχ)εχρ ibx AF[(B—byA] 

Sts (a? εχ! (x/a)exp—a\pl,a>O 
. af exp—x? ntexp—1f? 


i te: (—1y"-! cosnmx/c _ es 
ie, f(x) F,(B) = [Ὁ fesmin pds » 


F, (x) 3mf (B) 
anf (x7!) F (Ap) 
χ᾽ sinh ax cosech yx tan” *[tan(}roy~ ')tanhGay-* *p)) 
χ᾽ texp—ax,a > 0 tan~'a~*B 
exp—ax βία" + B?)~* 


fc) F(B) = §& f(xjeospxdx 
F(x) tn f(B) | 


exp—ax? Anta her β Μὰ ἘΣ 


ny lies Ὁ 


cosh ax sech yx cos nay +cosh my" 
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